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ABSTRACT 


A digital computer program to provide the control 
engineer with an interactive system for the analysis and 
design of linear continuous feedback control systems has 
been written in FORTRAN. The program, implemented on an 
IBM 1800 digital computer operating under the Time Sharing 
Executive system, required 9700 words of core storage plus 
disk storage amounting to about 400 sectors. 
The control systems design and analysis program 
(CSDAP) uses the block diagram, transfer function represen- 
tation technique, familiar to control engineers, for prob- 
lem formulation and subsequent data entry. The individual 
transfer functions are entered through the user's keyboard 
where they are immediately error checked before being 
included in the system definition data. After the trans- 
fer functions have been entered the user must then decide 
upon the technique or techniques to be used in the design 
and analysis of the control system and select accordingly. 
The user can specify one or more of the following options: 
a) The open or closed loop transient response of the con- 
trolled variable to a change in set point or load vari- 
able can be determined. Various performance criteria 
can be calculated from this response data. 

b) The frequency response of the open and closed loop can 
be calculated and the results can be displayed as a 


Bode, Nyquist or Log-Modulus (Nichol's) plot. In the 
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case of the Bode diagram the gain and phase margins 
and crossover frequencies are included. 

c) The root locus values are calculated for both posi- 
tive and negative values of the gain. The results 
are available as a root locus diagram or as a listing 
of the frequency coordinates and gain values. 

The calculated values are available as a display on 
the oscilloscope, as graphs drawn by the Calcomp Digital 
Plotter or as a tabular listing on the line printer or as 
any combination of these methods. 

The program has been designed with a strong emphasis 
on an interactive "conversational" mode of operation so a 
considerable amount of user assistance has been provided as 
part of the program. This assistance, which is available as 
a user selected option for all sections of CSDAP, minimizes 
the prior study required to utilize the program. The pro- 
gram allows the user to alter particular parameters of in- 
terest without re-entering the remainder of the parameters 
which remain unchanged. 

A complete Systems Manual for CSDAP is available so 
modifications to the existing program, for use on an IBM 
1800 operating under a different executive system or another 


machine, should not be difficult. 
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CHAPTER I 


INTRODUCTION 


The conventional use of a digital computer for ana- 
lysis and design problems suffers from two major disadvan- 
tages. 

First, the usual form of supplying information to the 
computer programs is by means of cards or tapes. Those 
which must be prepared according to rigid formats dictated 
by the program are often alien to the design engineer. 

Even in the case of programs in which an effort has been 
made to simplify the required inputs that are used, the 
designer is usually plagued by input data errors until he 
develops a certain familiarity with the particular programs. 
Errors due to improper input data preparation disrupt the 
design iteration cycle and frustrate the designer. 

The second and more important disadvantage is the 
delay, due to a long turn-around time, in getting the com- 
putational results back to the designer. This delay under- 
mines the effectiveness of the iterative design process 
and makes it impossible for the user to try out new ideas 
rapidly. 

A means of bypassing these man-computer communica- 
tion problems lies in the use of on-line computer consoles 
which have graphical input and/or output capability. 

Conversational programs using such equipment allow 


the user to input information into the system that may take 
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the form of command functions or data. The solution is 
then available in the form of a graph or table of print- 
ed results. Modifications to the problem parameters 
based on the results of the previous calculations can 
be entered and the problem rerun. 

It matters little to the user how efficiently 
such conversational flexibility is achieved, as long as 
the programs are reasonably responsive and easy to use. 
From the systems point of view, there will always be a 
serious concern for the amount of core memory and bulk 
storage required and the execution speed. 

This thesis describes how a set of programs that 
have been developed provides a versatile means of using a 
digital computer to analyze and design linear control svs- 
tems. The programs are written in the FORTRAN IV language 
and were implemented on an IBM 1800 computer located in 
the Department of Chemical and Petroleum Engineering. 

The IBM 1800 computer operates under the Time 
Sharing Executive (TSX) operating system which allows a 
background or nonprocess program to utilize any computing 
time thatis not being used to monitor and control on-line 
processes. Peripheral devices on the computer such as 
disk bulk storage, card reader/punch, line printer, key- 
board/typewriter, digital plotter, or storage display 
oscilloscope allow the analysis programs to operate under 


the TSX monitor in a conversational mode. 
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The total system is organized in much the same way 
that a control engineer might approach the problem. First 
the transfer functions of all the component blocks in the 
control loop are found and entered into the system. Once 
all the transfer functions are in storage, a grouping of 
all the transfer functions is made and used to form the 
open loop and closed loop transfer functions of the system. 
Subsequently, the time response, frequency response, and/ 
or root locus plots can be produced at the user's option. 
These topics are covered in more detail with examples in 
the chapters that follow. 

Although this thesis briefly describes most of the 
subroutines used, it does not contain extensive program 
documentation. The Systems Manual that is available can 
be consulted for the organizational flowcharts, program 
linkages and program listings. A User's Manual is included 


in Appendix A. 
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CHAPTER II 


LITERATURE SURVEY 


There are a number of methods that may be used to 
perform an analysis of the behavior of a control system. 
Although themajority are graphical techniques, they can 
be calculated analytically with the use of high speed 
digital computers. The following sections will discuss 
some programs and analytical methods that have been used 


previously. 


RatiControl Systems Analysis Programs 


The last decade has seen a number of programs writ- 
ten to handle the laborious task of control systems analy- 
sis. Cummins (12) developed a generalized procedure for 
synthesizing automatic control systems using some of the 
fundamental concepts. The procedure was programmed for 
the Ferranti-Mercury and the IBM 7090 computers. Computa- 
tions are performed after the system is described by a 
linearized set of equations. 

In the early sixties with a great deal of interest 
focused on the aerospace industry IBM (15) developed a 
program for its 7090 computer. This control systems ana- 
lysis program was written primarily to accomplish real- 
time flight simulation by digital computational techniques. 
The program provided a convenient method of obtaining a 


root locus for linear, continuous or sampled-data systems. 
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Input data can be in the s-transform and/or z-transform 
notation. Optional features included the conversion of 
a polynomial in Laplace notation to sampled-data or 2Z- 
transform notation and the computation of pole-zero loci 
and gain loci for a specified gain. 

Lofkrantz (24) in an investigation of this IBM 
program found that routines could be easily incorporated 
to allow the plotting of the Bode and Nyquist stability 
graphs. A method for calculating the transient response 
from the z-transform was also implemented. Some modifi- 
cations were also made to decrease the number of runs re- 
quired when calculating the closed loop transfer function 
and its subsequent crane eee response. 

The analysis and design of multiloop, multivariable, 
linear control systems requires the use of computers. 
Coffey (8) introduced algorithms for applying frequency 
response, pie ssoonscn and root locus methods Pomeuch 
systems with the aid of a digital computer. He describes 
a digital program AUTO, which uses the gradient optimiza- 
tion scheme to automatically synthesize compensation for 
multiloop continuous and digital control systems. This 
compensation is calculated to achieve a desired open-loop 
frequency response. The second part of his study is de- 
voted to a discussion of simulating complex systems by 
means of digital, analog, and hybrid computers with empha- 


sis being placed on the fundamental aspects of simulation. 
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A recent program developed by Agostinis (1) allows 
the control engineer to use a digital computer to assist 
in the analysis and design of linear, time-invariant con- 
trol systems or multivariable systems involving pure time 
delay elements. The program uses the matrix formulation 
of the state variable equation to produce the time domain 
solution for the output variable. The study gives a num- 
ber of worked examples, a user's manual and a good section 
on system documentation including detailed flow charts and 


program listings. 


B. Transient Response with Associated Performance Criteria 


The Laplace transformation is the most common means 
of solving the differential equations describing the 
performance of a control system. In order to find the 
transient response it is necessary to invert the Laplace 
transform. Corrington (10) showed that the transient res- 
ponse for a system described by the ratio of two polynomiais 
in the complex frequency domain could be calculated using 
a simple linear difference equation. Liou (23) formu- 
lates the corresponding state equation from the transfer 
function, derives a recursive formula for the solution, 
and gives the result for any desired interval of time. 

Chen and Yates (6) discuss a matrix formula for the inverse 
Laplace transformation. After the eigenvalues and coeffi- 
cients are substituted, and some simple matrix operations 


are performed, the inverse Laplace transformation of the 
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function being studied is obtained. 

The matrix formula by Chen and Yates has been im- 
plemented on a digital computer and some of the algorithms 
and the program listing appears in the book by Chen and 
Haas (4) and is also given by Chen and Shieh (5). The 
formula is general and is particularly advantageous for 
use on high-order transfer functions. Since many of the 
control systems studied in the process industries are of 
a moderately low order, the method of partial fractions 
could be efficiently used to invert the Laplace transforms. 
Ghtsim ofa tramicivaw.dtéchni ques (7 whieh, tht, 22.9 26, v28A10e29, 33) 
and with the help of routines (17) for polynomial manipula- 
tions and a Newton-Bairstow (9) approach to calculating the 
poles of the system, it could be implemented on a digital 
computer.. 

The pure time delay or transportation lag as it is 
called is frequently found in control systems. It cannot 
be expressed exactly as a ratio of polynomials. Consequent- 
ly some form of approximation is desirable before the tran- 
Sient response can be calculated. As some authors (11, 26, 
29) show it is frequently approximated by a function with 
a high order real pole far from the origin, or by a ration- 
al function of s that includes the first several terms in 
the power series expansion. This latter approach is referred 
to as the Pade approximation. 

Some compensation is usually necessary to improve the 


transient or steady state performance. The principle 
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difficulty is the establishment of a single criterion 
for performance that can be used for design. Several 
paperss (25,°27,°30)“and a number of. textbooks (13, 22, 
26, 29) have presented and compared several figures of 
merit for optimizing the system response. Several of 
these performance criteria have been included in the 


program. 


Ca 7hOO tt. LOCUS Analysis 


The graphical method of constructing the root 
locus of the closed loop poles, when the gain is varied 
from zero to infinity, has been well developed (14, 3l). 
An analytic approach for obtaining the root locus with 
positive and negative gain has also been described (3, 
32). These studies show how solutions for various sys- 
tems consisting of rational polynomial functions of s 
are obtained. A more recent study shows that the ana- 
lytical method merits utilization because of the flexi- 
bility inherent in the manipulation of algebraic expres- 
sions. In this study, Krishnan (21) demonstrates that 
arrays can be developed so algebraic expressions for the 
root locus and associated gain can be obtained from a 
knowledge of the open loop transfer function. Breakaway 
points on the real axis and the intersection with the 
imaginary axis can also be obtained directly from the 


arrays. 


i vets son Ree est 


a 


78 | | 


noitetixo eipnie 5°36" sisonie Retains gird TuSERVED | 
“Ef Ege RS: A@ERSD GOT Déau' Sc AMSG Se iditinddtiey oe 
Emer ery ‘@leeasuss 28" <aiwa"'s ae 
40 BstOpit isveves bots qtoS Bas besnéestg svsa “— 
86 £4698088 =) Senogher’ mMS3eys Sd Pate emisyoo wes oF 


7 “ 
bay 


eo tae 
aa3 


iS SAa°HE BSby tsa? need: oven’ BP+8s 19S GoAsmeei wg 


i t 
Fy Sytgas a #5 ce eps eee - or 0" severe re te rg eg Te, See 
t > ned ee Pace =e ae Y 2 yas 3 Be Week ao Uae pl Rae di Vek Lad viky) Wee Seer 
bg ’ , ‘ 
: : 
* 
ot a - : 


vg ‘ , ; . : . i ‘ = 
ty". ee Fry oy ee ee Pl Yee a a ms ty ee eos Pm Prat se. . py 
mod ye NY "2 SINS 7s ie!) Sey oh * weed "SG .7 G tae ‘ S| a 


eley AR es. 


vitnebartig ed Bition erorsoaat Le iPyeeeeeeees 
Pres qtQOk edt. paitouztene9,.i¢. bodten. Pspidaayp eet, 
oi ad D@kFeN, et aise eds eedye yesieq gqeol; peeps edt. AP ERM 
othdh 994). Degohoved.. fen deed.86d, eer © 
—— eyood 290% .9rt..padaiatda 2ol Apagugge Pf 
+£) Sadiunaad aeed gel 6, Red Fee oMaiepan: bos. OG 
Bal BUGLISY, 303. aqgisulos..wos. sWORG: ROSDUSS ORRET ats | 
Jon@> 20; Baodsagyd. igsdmoaylog, taapitsy: 2e.paddqtegon.t meg 
‘aay Re Oct pads BWORS yausa.. £45997, SORE: of i : : 
. thwelt..9A4.0..98ussed- not Si bate aideied alban 
ao “BetaKe ,2assdogts,.to, sto apie anbigih deen : 
i+ SHS  8ed69 enomeb.. (1S) aerate, seh rie | 
so9d8.298 anohocatgxs, obaxdapls fails ace R83 ays 
‘ni fO2? ReRbSt¢o..e4 geo akepyagseigoRes 2 
Mews Agen ...asigny3 ASRBASTT Good, ego wt 8 


eft dtiw noitosersjni stit panier sei: 


rh 
-_ 


7 Pont aa Si ae seca cal in i: 


7 ma 


rts Pie a bilby teak Bir 


D. Other References 


During the development of the program a number of 
manuals associated with the use and operation of the IBM 
1800 data acquisition and control system were consulted 
(15-20). 

Benunbear Of textbooks (2, il, bo, 227 .c0e.20).4awere 
utilized for methods and examples of the frequency response 


analysis techniques. 
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CHAPTER: LLL 


SYSTEM ORGANIZATION 


ATe Control System Representation 


One of the first steps when analyzing a physical 
system is to derive a mathematical model which describes 
the characteristics of the system. This representation 
of the interrelationships between system components allows 
the control engineer to use the available mathematical 
techniques to study the behavior of the system. 

The analysis of linear feedback control systems 
is facilitated by the use of transfer functions and block 
diagrams. A block diagram is usuaily described as a func- 
tional layout of an entire system in which the transfer 
functions of the components are shown in the blocks and 
the flow of signals is shown by lines connecting the 
blocks. Each block is an independent unit in that the 
connections between the blocks do not affect the transfer 
function inside it. Aithough the block diagram is used to 
provide a basis for an analog computer solution, it is 
also useful in the formulation of a problem which is to 
be solved numerically using a digital computer. 

The transfer function is obtained by writing the 
differential equation relating the output to the input of 
the block diagram, replacing the operator d/dt by the 
complex variable s, and solving for the ratio of the output 


to input. It is one of the most general and convenient 
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means of characterizing the dynamics of a physical system. 

The combination of the block diagram and the trans- 
fer function of a physical system provides a pictorial 
representation of the cause and effect relationship between 
the input and output of the system. 

Although all systems that are characterized by one 
input and one output may be denoted by a single block 
connected between the input and output, the advantage of 
the block diagram concept lies in the fact that feedback 
control systems are composed of many non-interacting ele- 
ments whose transfer functions are determined independently. 
Therefore, an entire system may be represented by the inter- 
connection of the blocks of individual elements, so that 
their contributions to the overallsystem performance may 
be evaluated individually. 

In the study of single-loop control systems, certain 
transfer functions are of special importance. In the 
general case shown in Figure 3-1,G(s) is referred to as 
the "forward" transfer function while H(s) is called the 
"feedback path" transfer function. The relation between 
the output, C(s), and the input, R(s), of this system can 


be obtained as shown in Equations 3-1 and 3-2. 


G(s) R(s) - G(s) H(s) C(s) (3-2) 


Cs) 


C(s) G(s) a 
Ris) -*1 + Gls) H(s) (3-2) 
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C (s) 


NOMENCLATURE 
R(s) - Setpoint or reference input 
C(s) - Output or controlled variable 
G(s) - Forward path transfer function 
H(s) - Feedback path transfer function 

G(s) H(s) - Open loop transfer function 
G(s) 
1 +6(s) H(s) - Closed loop transfer function 
1 + G(s) H(s) = 0 - Characteristic equation 


FIGURE 3-1 STANDARD FEEDBACK CONTROL SYSTEM BLOCK DIAGRAM 
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Since the transfer function in Equation 3-2 relates 
the output to the input with the,loop closed, it is refer- 
red to as the “closed loop") transfer function of the system. 
The relationship between the output and the input of the 
system through the forward path transfer function when 
feedback action is removed by opening the path between H(s) 
and the comparator given by the transfer function, G(s) 
H(s), is known as the "open loop" transfer function. The 
denominator of Equation 3-2 is referred to as the "charac- 
teristic equation". 

Seldom is it possible to represent a control system 
with the simplified block diagram of Figure 3-l. The con- 
trol engineer is normaily interested in analyzing the be- 
havior of a system for disturbances in load variables and 
the effect of varying control parameters on the overall 
system performance. In reviewing the complexity of process 
control systems that might be studied, the author finally 
decided on the general system representation shown in the 
block diagram of Figure 3-2. 

Once a control system can be represented in the same 
form as this block diagram, the next step is to determine 
the "overall" transfer function that relates the output 
to the setpoint or to either one of the load variables. 

A symmetric approach for obtaining the overall transfer 
function for a setpoint change or load change will be des- 


cribed in Section C of this Chapter. 
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B. Master Program Description 


The following three sub-sections will describe 
the overall system organization of the programs and the 
manner in which they were implemented. The first sec- 
tion will be a discussion of why a conversational system 
was used while the other sections will outline how the 
executive controls the entry into different phases of the 
total system and how communication between these phases 
is implemented. 

1. Conversational Operating System 

Conversational programs were prepared to use ter- 
minal input and output facilities for carrying on a dialogue 
with users. User programs will always require some form 
of input from the user; either command functions or data. 
This is an obvious source of error because of lack of ex- 
perience, typographical mistakes, or transmission errors. 
It is most important to provide assistance to the user in 
guiding proper input, reporting detectable errors, and 
enabling corrections. This later facility is usually re- 
flected by immediate syntax checking. 

However, there is one difficulty with providing 
such input help to a user, namely that depending on the 
expertise of the user, it may or may not be appreciated. 
If a user is very adept, excessive verbosity on the part 
of the program will be considered a severe annoyance and 
will tend to make the operation slow. If the user is a 


novice or highly subject to making errors, it is appropriate 
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16 
to insure proper input. 

1.1 Optional User Assistance 

By recognizing that a user may optionally need 
comprehensive input help, a conversational system can 
satisfy the real-world spectrum of users. This fea- 
ture has been provided in the present version of the 
program by means of data switches. 

When operating the executive program checks one 
of the sense switches (determined at system loading time). 
If this switch is in the "OFF" position, the IBM 1800 
console data switches will be used. If the switch is in 
the "ON" position, the set of digital input switches will 
be used. 

These data switches when set to the "ON" position 
cause the program to list the possible alternatives or 
options for the requested function parameter. Table 3-1 
lists each section of the analysis program that individ- 
ual switches control. As the user becomes familiar with 
various sections of the total program, some of the switches 
would be left in the "OFF" position. 

1.2 Recognition of Input Errors 

Error processing is one of the most important fea- 
tures of a conversational program. Those errors or prob- 
lems occurring during execution which require some decision 
making or correction can be resolved because the user is 
available at the terminal. Programs operating in a batch- 


run environment usually assume that a serious error will 
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TABLE 3-1 


MEANING OF THE DATA SWITCH OPTIONS 


Program Section Controlled 
LIST CONTROL DIGITS AND SWITCH DATA 


LIST TRANSFER FUNCTION DATA AFTER INPUT 


INPUT-OUTPUT VARIABLE AND FORCING FUNCTION 
SELECTION 


LIST PARTIAL FRACTION DATA 

REPLOT, AND PLOT OPTIONS FOR TIME RESPONSE 
LIST TIME RESPONSE DATA ON THE PRINTER 
PERFORMANCE DATA OUTPUT DEVICE OPTIONS 
FREQUENCY RESPONSE DATA OPTIONS 

BODE PLOT OPTIONS 

LIST FREQUENCY RESPONSE DATA ON THE PRINTER 
NYQUIST PLOT OPTIONS 

LOG-MODULUS PLOT OPTIONS 

ROOT LOCUS PLOT OPTIONS 


LIST ROOT LOCUS DATA ON THE PRINTER 
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cause a job abort because the user is usually not pre- 
sent. Even if the user was present he would not be 
allowed to intervene because it would hold up the job 
queue excessively. 

In the control systems analysis program being des- 
cribed here the data is error checked immediately after 
the user completes the input operation. If there is an 
error, an appropriate message is typed out in red print 
and the user is requested to re-enter the data. 

A typical example would be the case where an error 
had been made by the user when specifying the end time and 
subsequently the number of subdivisions for the time axis. 
The correspondence might have appeared as follows with the 


user response underlined. 


ENTER TOTAL TIME 


LK 
“Al ERROR IN INPUT - TRY AGAIN 
25 
"a ENTER NUMBER OF TIME DIVISIONS, MAX. = 500 
600 
DATA WAS NOrCONSISTENT, TRY AGAIN 
60 


The computer subroutine first found that the user 
had entered non-numerical characters to specify the end 
time and so requested the data again. The number of time 
divisions after it is entered is first checked to make 
sure it is a positive integer and then it is compared 
against the limits and in the above case it was greater 


than the maximum allowable. 
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19 
1.3 Convenience in Inputting Data 


Although the majority of the programs were written 
in FORTRAN, the author found it necessary to write a num- 
ber of subroutines associated with the input of data to 
the program in IBM 1800 assembly language (19). 

The first of these was to enable the user to enter 
a single numeric digit. from the typewriter keyboard. This 
was necessitated because many of the options or branch 
decisions could be coded using the numeric digits 0, 1, 2, 
--e, 9. The conventional method meant that the desired 
numerical character must be typed followed by the end of 
field character. This extra character provided for addi- 
tional delay or errors to occur. 

The calling sequence of this subroutine INCHR which 
is listed and flow charted in the systems manual is given 
below. 

CALL INCHR (IDEN, INMBR, ILOW, IHIGH) 

It allows a numerical character INMBR to be entered through 
the keyboard on logical unit number IDEV. The other two 
parameters ILOW and IHIGH provide an allowable range for the 
inputted character. The subroutine reads a single character 
from the keyboard, checks that it is a numerical character 
and if it is that it does not lie outside the range pro- 
vided by the limits. If either of the above errors occurs 
the routine types a message and requests it again. 

Another subroutine FFINP allows the user a less 


restrictive method data input to his program than is 
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available in FORTRAN. The subroutine allows integer 

or real variables to be read in as single values or as 

a vector of values. The variables are entered separated 
by a space or comma delimiters. The input field can con- 
tain any number of leading and/or trailing blanks but 
cannot contain any embedded blanks within the variables. 
Alphanumerical data can also be read into the integer or 
real variables under the A2 or A4 format respectively. 

In the case of an alphanumeric data entry the field begins 
at the first non-blank character and continues until the 
field is filled or a common delimiter is encountered. A 
space is not a legal delimiter for this data entry. The 
subroutine is fullydocumented in the Systems Manual. 

The final subroutine written in assembler language 
was TWAIT. This routine when it is called, loops inter- 
nally until the device specified in the calling sequence 
is not busy. It is required because a keyboard entry is 
of higher priority than outputting a message and a request 
for input might be made before the description of the 


alternatives was listed for the user. 
2. Executive Program Organization 


A conversational program permits the user to change 
his mind by Stopping the execution or starting the pro- 
gram over again, and even returning along the calculation 
path to a previous logical breakpoint. Such flexibility 


Supports the dynamic on-line user environment. 
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The Control Systems Design and Analysis Program 
has thus been divided into ten logically self-contained 
sections as shown in the simplified program flowchart 
Figure 3.3. This program requests the user to enter the 
CONTROL DIGIT for the next section of the program to be 
executed (see Table 3-2). The list of these options can 
be obtained by the user at the start of a run and option- 
ally during the program execution by having data switch l 
in the "ON" position and entering a CONTROL DIGIT of 0. 

The CSDAP System is designed so that as each sec- 
tion is completed it returns to the executive. This means 
the calculational path can be changed at any logical break- 
point which is the successful completion of one of the sec- 
tions. 

For example we will assume that the user has com- 
pleted entering the necessary transfer functions, calcula- 
ted the inverse Laplace transformation, and displayed the 
time response on the storage scope. It is possible to 
have the same results plotted onto the digital plotter by 
repeating only the programs associated with CONTROL DIGIT 
number 3. Because the plot is just a replot of the data 
previously calculated for the scope display it is only 
necessary for the user to enter the plot device, plot 


direction, and the title for this graph. 


3. Communication Between the Program Sections 


The communication between the programs and subrou- 


tines was achieved by using FORTRAN COMMON and disk FILES. 
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FIGURE 3-3 SIMPLIFIED FLOW DIAGRAM OF THE PROGRAM 
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TABLE 3-2 


LISTING OF CONTROL DIGIT FUNCTIONS 


Control 

Digit Program Section Entered 
0 LIST CONTROL DIGIT AND DATA SWITCH OPTIONS 
uy INPUT TRANSFER FUNCTION DATA 
2 CALCULATE INVERSE LAPLACE TRANSFORM 
-3 PLOT THE TIME RESPONSE DATA 
4 CALCULATE THE PERFORMANCE CRITERIA 
5 FREQUENCY RESPONSE CALCULATIONS 
6 ROOT LOCUS PLOT 
7 LIST TRANSFER FUNCTION DATA 
8 INITIALIZE TRANSFER FUNCTION AREA 


9 CALL EXIT 
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Their use will be briefly described here with full docu- 
mentation in the Systems Manual. 

The COMMON area is defined as four blocks and has 
therefore kept the COMMON statement in each program as 
simple as possible. The first block ISTOR is an array 
of eighty integers that contains variables frequently 
used by all of the program sections such as the logical 
unit numbers of the input/output devices, file numbers, 
array pointers, data switch settings and the names of 
the transfer functions. It is initialized once at the 
beginning of a run by the mainline CSDAP for any variables 
that remain constant during the run such as the logical 
unit and file numbers. If a variable is required ina 
program it is available by equivalencing an integer to 
the correct element of the ISTOR array. 

The other three blocks of COMMON are used to store 
the transfer function data and appear in dimensional 
form in the program as IC (2, 40), ICONT (200), and DATA 
(400). Their use will be discussed further in Section C 
of this Chapter. 

A few disk files are used in the system. The main 
use of these files is for saving the COMMON area and for 
storing calculated data for replotting the time response, 
the frequency response diagrams or the root locus diagram. 
Because the size of core available to the nonprocess user 
is under ten thousand words, it was found necessary to 


overlay the common area when not directly working with the 
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transfer function data. This area was saved and restored 


into a file on disk by the use of subroutine DATSV. 


Cc... Transfer .Functions 

1. Block Diagram Transfer Functions 

The block diagram that appears in Figure 3.2 has 
been used as the basis for the transfer function entry 
into the program. The control system is composed of two 
control loops and nine blocks with the reference or set- 
point and two load variables as inputs and one output 
variable. 

To facilitate easy block diagram manipulation, the 
transfer function of each block is entered as a ratio of 
polynomial terms in the complex variable s. A general 


form taken by such a transfer function would be 


gq a foul Ue 
s + a,) (ags + ass eters ire eta) 


p Ret 
+ bis + b,) (bes + bs ets cease) 


a 
(ajs + a) (as + ay 


P3(s) = 5 


(bjs + b5) (b3s 


(353) 


The transfer function (Equation 3-3) is identified 
in the program by the two alphanumeric characters P3 which 
represents process block number three. A similar designa- 
tion is used to represent the other individual block trans- 
fer functions shown in Figure 3.2. 

2. Transportation Delay 

The pure time delay or transportation lag as it is 


sometimes called is frequently found in control systems. It 
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cannot be expressed exactly as a ratio of polynomials 
as was the case for the other transfer functions. It is 
therefore necessary to make some form of approximation 
which preserves the all pass nature of the pure time delay 
and gives an essentially linear curve of phase shift ver- 
sus frequency for the range of frequencies being investi- 
gated. 

There are several means of achieving the required 
transfer function for the time delay. One method of approxi- 


Ts 


mation involves the replacement of e by a rational func- 


tion of s such that the first terms in the Taylor power 


series expansion of ane appear in the denominator (Equation 
3-4) 
1 
e TS = to = Sr a ey a i Te 3 (3-4) 
e Lies FS / 20 B/S. OE ere 


This approximation is reasonably accurate for values of 
wt<l but for wt>l both the magnitude and phase shift differ 
sharply from the exact representation as shown in the Bode 
plot’ (Figure 3-7). | 

The Pade approximation (26) for the time delay is 
used to match any required number of coefficients in the 
power series expansion, by a ratio of zeros and poles. The 
more accurate approximations require more poles and zeros. 


A general formula for the Pade approximation appears in 


Coughanowrand Koppel (11). 


as : .” Vie 


7 


afeimonyiog Jo ortst 6 26 tala a eiutaniexe od ore 


- - 


to? sesD snd swe 


”-~ 
4 
~~ 
© 
] 
4 
|e 
o) 
ti] 
-_ 
~_ 
=e 


noigsmixorqas to mrot smo2a sasm od Yesses0en sxotexeds 


ysieb amis stugq sft to siuz6n eesq Ife sm asvisee%g Aokw 
_ 


“sav tiinde seera to svavo azsonil yi lsijnesase as eS vip ‘bas 
ee 
-tt2zevoi poised satonsupssi to spnsi sdt xo yoneupes® ave 


o_ 7 
Seiiups1 edz palveiros ito enssm [stevse ous Szen? : 

: aft . 
qqe to botitem enO .yeheh emts st sot nolssaut =e8 
‘ CP 2 
“onui Lanoivs1 s vi 9s to tnenmecsiaqst ait eevioynt nose 
pe . 
‘ ; 
xewod t0olvaT eit at-erties sett? eft 2add dove «a. to mols 


i 
¥ 


; . . ; i 
nolssupS) yossnisenaS eis al tasgas a xo MOLENSGXO eo!x08 


to aesuisv 102 ststuoNS videror ‘BST Bl: noivemixozggs 

steTtTLb tiide sasdq bas ebuctinpem ods dtod f<su 2102. stud J 
gbod sit mi nwode es molistsinsesiqe1 tosxs eft mozt yiqus 

(THE suvpla) se - 

et yAlob emit srit tot (08) noisemixevgqs e669 ont 


sat at etnstortisos to us 


“t+ 


: * 
mun bottupst yns dosem oF Bese 


:% 


—— 


siT .aslog bus sotes to oizs1r s yd \tokeasaxe eolzes 
* al 
8018s bas eelog stom srivpst snotssméxoxags oaswions | eit0m 
a i” 
bral ‘easeqqs noltemixo1ggs obsd ents 392 tumor 1 fsx — aoe r 
es 


yrs 


7 ; i 


aa lie 


o 
aiid 


7 WY 


palin 


ath 


The first four orders of the Pade approximation 


are listed below: 


—TS 1 = +s/2 


e = 1 eae AD (Firet Onder) (3-5) 
Bet ite/2 +, eB / ie 
e = Zz) (Second Order) (3-6) 
Poets /O + TS / 12 
2, 1.-.18/2°+t7s*/12 - 17s°/120 
e = a) 35 (Third Order) (3-7) 
te as/2-+t_ ss fl? + ss" /120 
= rue dice 18/2 + 4°S7/12,2) 12°8°/120» 4,.1/8°/1680 
as = a) 5 aed (Fourth Order) 
TS (Oe to ee Al) eee 120. ten 8 ./168,0 


(3—8) 


For these approximations the magnitude of the frequency 
response agrees with the true magnitude exactly for all _, 
and the phase shift is a better approximation over a wider 
frequency range than the power series approximation (Figure 
3-6). 

The time delay approximations are handled by a sub- 
Poul ne DLYTM that is entered from the transfer function in- 
put section. Instead of a legal transfer function name as 
shown in the block diagram (Figure 3-2), the letters TD are 
entered for deadtime. The DLYTM subroutine is then entered 
and requests the user to enter further information describ- 
ing the time delay approximation. The time delay can be 
placed in the forward path-block P3 or in the feedback path- 
block Hl. Once the rational polynomial approximation is 


calculated it is multiplied by the existing contents of the 
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respective block P3 or Hl. This makes it possible to 
have another transfer function already present for this 
block and dictates that the time delay must be the last 
transfer function entry for either the P3 or Hl blocks. 

The two methods previously discussed for approxi- 
mating the time delay are both available in the program. 
For the power series method up to a fourth-degree Taylor 
series expansion can be used in the denominator. The 
Pade approximation can take any of the forms described 
by Equations 3-5 to 3-8 inclusive. Figure 3-4 shows the 
time response of both the first and second order Pade 
approximations to a unit step forcing function. This 
figure compares well with a similar figure given by 
Buckley (2). 

The time response for the fourth-degree power 
series and Pade approximations is shown in Figure 3-5 
where each of the time delay transfer functions had a unit 
step forcing function applied. The Bode plots of both 
these approximations is given in Figures 3-6 and 3-7. 

.3. Transfer Function Data Entry 

Thetinput of trarsfer~function data’ can best be 
explained by giving an example with a discussion of the 
topic. The coefficients for each polynomial of the trans- 
fer function (Equation 3-3) are contained within brackets 
and appear ordered from the smallest to the largest power 
of s. The distinction between numerator and denominator 


polynomials is made by the insertion of a division sign 
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FIGURE 3-4 Pade Approximations for Time Delay 
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after all numerator polynomials have been entered. 
The input for the transfer function given in 


Equation 3-3 would take the form: 


P3 (45,8) (a57a,7A43) (.--)/(b5,b,) (be rb, ,b3) (..-) (3-9) 


This data is read in by a subroutine INPUT which per- 
forms a legality check on the first two alphanumeric 
characters to insure that they represent one of the nine 
blocks or 'TD' for time delay or transportation delay. 

If the data passes this check, it is nen passed 
to a subroutine GDFIB which will analyze the rest of it 
for the numerator and denominator polynomials. The data 
is then transferred by subroutine MDATA to the transfer 
function storage area and substituted for any earlier 
version of the P3 transfer function. 

An error message is typed to the user if the enter- 
ed transfer function name is not legal or if an error is 
found in the coefficient data upon further processing by 
GDFIB. For the transfer function name error the INPUT 
subroutine exits back to the executive program which then 
requests the next control digit to be entered. This there- 
fore provides a convenient method for terminating the in- 
put section of the program. The pushing of the "end-of- 
field" character is interpreted as an illegal transfer func- 
tions name and is used for ending the inputting of transfer 


function data. 
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In contrast to the illegal name error an error in 
the input data is handled differently. Because the user 
has successfully entered a legal transfer function name 
the program assumes that he really wants to enter this 
transfer function. Therefore the program after typing 
an error message, requests that the data be entered again. 

To configurate a typical control system the trans- 
fer functions of a number of the blocks must be entered 
into the system. The order in which these transfer func- 
tions is entered is irrelevant. An important point to 
keep in mind is that whenever a transfer function is 
entered correctly it replaces the earlier version of that 
transfer function. One exception to this was noted in this 
Chapter (Section C-2) where it was explained that the time 
delay approximation was multiplied by the transfer function 
previously entered for me P3or Hl blocks of, Figure 3-2. 

4.- Storage of. Transfer Function Data 

The transfer function data takes a form in storage 
Similar to the input form shown in Equation 3-9. The poly- 
nomial data is stored as an array of coefficients ordered 
from the lowest to the highest power of s. Since the size 
of each numerator and denominator polynomial is not con- 
stant some form of data storage must be used so that the 
area used by each transfer function can float depending 
upon its own size. This means that a storage area equal 
to the largest transfer function will not be required for 


each transfer function of the system. 
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A method has been utilized that keeps the coeffi- 
cients of all transfer function polynomials in a contigu- 
ous string in storage. It makes use of three of the blocks 
of COMMON IC(2, 40), ICONT (200), and DATA (400), as men- 
tioned in Section B-3 of this Chapter. The IC array is 
used to store pointers for the other two arrays. The 
ICONT array stores the number of polynomials followed by 
the degree of each polynomial term that appears in the 
numerator or denominator of each of the transfer functions 
involved in the system. Similarly the DATA array contains 
the coefficients associated with each of these polynomial 
terms. 

The system handied a total of twenty transfer 
functions composed of the nine block transfer functions of 
Figure 3-2 which the user enters and the rest such as the 
open and closed loop, the forward path, and the overall 
transfer functions being calculated by the program. There 
will therefore be forty numerator and denominator storage 
areas required. 

The program associates a number ITF between one and 
twenty with each of the twenty transfer functions in- 
stead of using an alphanumeric representation similar to 
the user.” The-values—of-i¢—CL;7—2" 11TF-1)and-IC (1,2* ITF) 
give the first entries in che ICONT array associated with 
the numerator and denominator respectively, of the ITF 
transter function. The values of IC’ (2,2*ITF-1) and IC 
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vector for the coefficients associated with the numerator 
and denominator polynomial terms respectively. 


The method discussed above can be illustrated in 


the following example. 


Block Transfer Function: 
P2(s) = a 
Data Entry to. Program: 

B2e(6025 bull Onrhb/ (nel) b2, tba (Sfel) 


Storage in Arrays: 


Assume that 1) P2 is the fourth transfer function 


2) iC Cig ek 


3)) SECa 257). = 40 
Then Ic (1,8) = 18 
EC» (2)7-8)s = 43 


TCONT array ~-starting at FCONT (15) “contains 
Ze er eo pe PL 
DATA array starting at DATA(40) contains 


zl Ler Tyo, UU, ere OyeseOye ies, O07 130 


5.. Listing of Transfer Function Data 


Each transfer function can be listed automatically in 
the abbreviated form it will take in the storage array. This 
listing will be given if it has been entered correctly and, 


if the correct data switches have been set previously. 
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Another listing of this data can be obtained at any 
time by the entry of the CONTROL DIGIT 7. This option is 
used for checking previously entered transfer functions or 
for getting a complete listing of all the transfer func- 
tions for documentation purposes. The listing can be made 
on either the user's typewriter or the line printer by 
the subroutines LISTD and LSSTF. When a complete listing 
is made onto the printer, the program requests a title to 
be entered that will be associated with the listing. 

The various forms that the transfer function list- 
ing can take will be illustrated in the Appendices. 

Ors Manipulation of Transfer Function Data 

6.1 Initialization of the Transfer Function Data 

Arrays 

The first time the program is executed or upon 
the user's command (CONTROL DIGIT = 8) an initialization 
function is performed on the transfer function area. A 
subroutine INTFA is executed which sets the order of numer- 
ator and denominator polynomial to zero. Each coefficient 
is made equal to a constant of 1.0. This means that all 
transfer functions would have a value of 1.0/1.0. 

6.2 Inner Control Loop 

There is one exception to this procedure. It 
is the case for the inner control loop whose numerator is 
initialized to a constant value of 0.0. This value, 
which means that there is no feedback for the inner con- 
trol loop, remains until a different transfer function is 


entered for H2(s). 
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6.3 Open and Closed Loop Transfer Functions 


The open loop and closed loop transfer func- 
tions of the overall system can be calculated once all 
the polynomial terms in the numerators and denominators 
of the block diagram transfer functions have been entered 
into storage. The routine OLOOP is executed first and 
performs a block diagram manipulation in a manner similar 
to that normally employed in block diagram simplification. 
The inner control loop is first checked for feedback. If 
it has feedback, it is converted to a single block in the 
forward path by the subroutine CANON. The transfer func- 
tion calculated for the inner control loop is then multi- 
plied by the other controller and process transfer function 
blocks to obtain the forward path transfer function which 
is then transferred to the storage area. The forward path 
transfer function is multiplied by the feedback path trans- 
fer function H1l(s) to complete the calculation of the open 
loop transfer function which is also saved in storage. 

The closed loop transfer function is calculated by 
the subroutine CLOOP. Previous to entering CLOOP another 
routine CLTYP is executed which requests the user to enter 
the type of response that he desires because the closed 
loop. transfer function will be used to calculate the tran- 
sient response. The program allows the user to substitute 
the open loop for the closed loop transfer function or to 
calculate the closed loop transfer function as the ratio 


of C(s)/R(s) or B(s)/R(s) as defined in Figure 3-2. This 
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option was included because although it is common to 
refer to C(s) as the output variable, this variable 
usually cannot be measured and the control signal is that 
from the measuring element transfer function Hl, namely 
(e429 Be 

If the closed loop transfer function B(s)/R(s) is 
desired, it is calculated by putting the numerator equal 
to the numerator of the open loop transfer function and 
the denominator equal to the sum of. the numerator and 
denominator of the open loop transfer function. To cal- 
culate the closed loop transfer function C(s)/R(s) a 
slightly more complicated method must be used to manipu- 
late the forward path and open ioop transfer functions. 

Extensive use of the polynomial manipulation rou- 
tines of the Scientific Subroutine Package (17) was made 
in calculating the open and closed loop transfer func- 
tions. Two other subroutines were written for transfer- 


ring transfer function data. GETTF determines a single 


transfer function for a particular block while TFMVE trans- 


fer a transfer function back to permanent storage. 
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CHAPTER IV 


TRANSIENT RESPONSE CALCULATIONS 


A. Laplace Transforms 


Application of Laplace transforms to linear dif- 
ferential equations with constant coefficients yields 
an algebraic equation in another variable s. The solu- 
tion of the differential equation is effected by simple 
algebraic manipulations in the s-domain. To obtain the 
desired time solution, it is necessary to invert the 
transform of the solution from the s domain back to the 
time domain. The Laplace transform method is used exten- 
Sively and the pertinent theory can be found in many books 
erm eedback;, control theory’ (11, 13; 22, 26, 29). 

Two important advantages of the Laplace transform 

method are: 

1) initial conditions are automatically incorpora- 
ted into the solution for any arbitrary forcing 
function, and 

2) solutions to problems involving simultaneous 
linear differential equations with constant co- 
efficients can easily be obtained. 

The usual practice in the analysis of control systems 

is to take the Laplace transformation of the original equa- 
tions provided the initial conditions for all variables are 


zero. If the variables in the original equations do not 
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have zero initial conditions a change of variables is 
made before taking the Laplace transformation. The 
transformed equations are then combined and re-arranged 
and finally the inverse transformation of the final ex- 
pression yields the desired solution for the transient 


response. 


B. Inverse Laplace Transforms 


The process of inverse transformation constitutes 
the major portion of the program for calculating the 
transient response. Even so, the procedure is orderly 
without using trial and error solutions, and the result 
yields both the steady-state, and transient solutions for 
given initial conditions. 

The inverse Laplace transform may be evaluated by 
means of the inversion integral and the calculus of resi- 
dues, but a much more expedient method is to use the par- 
tial fraction method. 

ite Eaucian Fracrion Metnod of inverting Laplace 

Transforms 

For simple transfer functions the inverse trans- 
forms areobtained by inspection. The more complicated 
transformed expressions require that a table of Laplace 
transforms be consulted to obtain the inverse transform 
function. If the expression whose inverse transform is 
desired is not in the particular table available, the com- 


plicated expression must be simplified into a sum of simpler 
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expressions, and the partial fraction method applied to 
perform the inversion. 

The denominator of the overall algebraic trans- 
fer function F(s) must first be factored to yield an 
expression of the form shown in Equation 4-1. 


eS). N(s) 
tapos Dis) (s+ $1) (s + 55) CET CE Ss.) (4-1) 


For an nth order denominator polynomial Dis), there must 
be n factors or roots. The values are the poles of F(s) 
or the values of s in the complex plane for which F(s) 
becomes infinite. These poles may be real, complex, or 
imaginary. 

If all the poles are different, the partial frac- 
tion expansion takes the form 
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where Ki,  =Atye2, ...'h ate constants. Once the K,'s 
are evaluated, the inverse Laplace transform may be taken 


term by term with the result being 


ree nek: a-Snt 


F(t) = K + K, + 6. + KL (4-3) 


1 


Subsequent sections outline how to perform the partial 
fraction expansion for different pole combinations using 


a digital computer. 
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2. Calculation of Partial Fraction Coefficients 
The response transfer function F(s) can be expres- 
sed as the ratio of two polynomials P(s) and Q(s). In this 
program, these polynomials must satisfy the following con- 
ditions: 
i) the degree of P(s) must be lower than the degree 
of Q(s), 
ii) these two polynomials must not have a factor in 
common, 
iii) the coefficients of the two polynomials must be 
real numbers, 
and iv) the coefficient of the highest power of s in the 
denominator must be equal to unity. 
The program assumes the user is aware of these conditions and 
so does not explicitly check for these conditions. 
The first step in performing the inversion by the 
partial fraction method is to factor Q(s) into linear and 
quadratic factors with real coefficients: 


Pts) P(s) 


F(s) = Q(s) - (s - S,)(s - So) ee S,) (4-5) 


These values s are called the poles of F(s) 


11 Sor cee Sy 
but are also sometimes referred to as the zeros of the 
denominator. Now having determined the K poles of F(s) 
it is possible to express F(s) as a series of fractions, 


where the number of fractions is equal to k. If there 


are no repeated terms, the function F(s) can be expanded as 


F(s) = == = = + = t+ ee s- 8, (4-5) 
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It is then necessary to evaluate the A constants and evalu- 
ate similar constants for repeated factors and complex 
factars of O(s)% 

The Heaviside Expansion Theorems as given in refer- 
ences (7, 15, 33) are of great utility in, systematizing 
the procedure of using partial fractions to invert the 
Laplace transforms. They outline how the coefficients 
for Equation 4-5 can be found for any pole of F(s) where 
it is real or part of a conjugate pair appearing in single 
or multiple order. There are four classes of problems en- 
countered depending on the zeros of the denominator Q(s). 

1) F(s) has first-order real poles 

2) F(s) has repeated first-order poles 

3) F(s) has a pair of complex conjugate poles 
(a quadratic factor in the denominator) 

4) F(s) has repeated pairs of complex conjugate 
poles (a repeated quadratic factor in the 
denominator) 

Each of the above cases will be discussed in some 
detail. — In general, the Laplace transform F(s) will con- 
tain combinations of these cases. Many textbooks (7, ll, 
13, 31) will give a more expanded description of the par- 
tial fraction method than will appear here. 

For the first case Q(s) has a linear factor s - a 
that is not repeated. If ¢(s) denotes the function which 
may be a quotient of polynomials that is left after re- 
moving the factor s - a from the denominator of F(s); such 
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F(s) = P(s) be o(s) (4-6) 


then according to the theory of partial fractions there 


exists a constant C such that 


'o- 

- 
6 

— 


+ H(s) (4-7) 
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! 
0) 
1) 


where H(s) represents the sum of the partial fractions 
that correspond to other linear and quadratic factors 
of Q(s) that appear once or are repeated. To evaluate 
the constant C, which is called the residue of F(s) at 
the poles = a, multiply both sides of the above equation 
by theterm (s - a) and as s > a it is seen that C = (a). 
The inverse transform of the fraction C/(s - a) is given 
by ¢ (a) sang The function ¢(a) can be evaluated from 
Equation 4-6 but 9(a) can also be calculated from the 
fact that 9(a) = P(a)/g'(a). This is an important fea- 
ture because all the constant C's corresponding to other 
linear non-repeated roots can be conveniently calculated 
from two polynomial evaluations for the numerator P(s) and 
the first derivative of the denominator Q'(s). 

Considering the second case where the polynomial 
Q(s) contains a linear factor (s - a) which is repeated 


n times. Following the same procedure as above it is 


possible to write 


» Bish = (s) 4 
Ps ie P45 is ee BA aL (4-8) 
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where ¢(s) is the quotient of polynomials obtained by re- 
moving the factor (s - a from the denominator of F(s). 
There are now more terms necessary in the partial 


fraction representation of F(s). The form now taken will 


be as follows: 


¢(s) A A AL 


= = + —————, + 1... + ——+ + H(s) 
(s - a) Sn s-a (s - ave (s - ayes 


G4 =9)) 


where the constants A. are not explicit functions of s and 
H(s) has the same meaning as in Equation 4-7. Again mul- 
tiplying both sides by the repeated term results in 


o(s) = A, (s ea) art A, (s - 5) BA weakens ena (s - ay ntl H(s) 


(4-10) 


where it is seen that hy 22 ¢(a) as s +> a. To find the re- 
maining n coefficients both sides of the above equation are 
differentiated n - r times with respect to s to isolate 

the Sona ants AL: When s > a in the resulting equation, the 


relationship for the A'S is found to be 


> (n * 1) (a) 
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Ay A rrett (n - xr) 


£ Er re ge eda), tren. 
(s - a x! in or pis! 
(4-12) 


Applying this relationship repeatedly as r goes from 
0 + n yields a summation of terms that will represent 
the inverse transform of ¢(s)/s - aynoe. 

In the following section imaginary linear factors 
will be discussed. Since it was stated at the outset 
that the polynomials P(s) and Q(s) would have real coeffi- 
cients, it follows that any imaginary linear factors of 
Q(s) resulted from a quadratic factor that could not be 
factored into two linear terms. The quadratic has the 
form 37 + as + 8, where aa - 48 < 0 and therefore results 
in two complex conjugate factors. By completing the 
Square, the original quadratic factor can be converted 
into a more useful form such as: (s - a)? + be where a 
and b are real numbers and s = a + ib are the complex 
factors of the original quadratic s? + as + § also. 

If the quadratic factor appears only once then 
the following relationship can be obtained 


R.( Sila dened & o (s) = As. + B 


F(s) = Stay Eee ot (8) 
Qis (s°= al? Ets. (g = a)? + 5° 


(4-13) 


where A and B are real constants and H(s) has the same 


meaning as in Equation 4-7. Multiplying both sides of 
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this equation by the quadratic factor and letting s > a 


+ ib gives 
¢{a + ib) = (a + ib) A+B (4-14) 


Uibe oy and d5 are said to represent the real and imaginary 
components of the complex number ¢(a+ ib) it can be seen 


that after suitable substitutions 


(s - a) dg, t+ bo 


5 + H(s) (4-15) 


1 
a(S We Be 


which has an inverse of the form 


F(t) =~ e** (9, COS bt + 6, SIN bt) + H(t) (4-16) 
Simplification of this expression can be achieved by con- 


verting to polar coordinates. It then would appear as 


at 


F(t) =~ |o(a + ib)| e** SIN (bt + 6) (any) 


> 
where © = argument of ¢(a+ ib) = ARCTANGENT (=*) 
1 
: 4 2 2 
g(a + ib)| = /o,° + 45 


A similar procedure to that for repeated linear 
factors must be followed in handling repeated quadratic 
factors. As the majority of the problems in the analysis 
of control systems can be handled by a quadratic factor 


that is repeated once, only this case will be considered. 


‘ 
' Bs 


+ 8 pnistel base rosos? oissetbsup, sit vel sobsaupe 
sevir. . 
(al=—5 ag +A (dt + 8) * (dk + BO 


yusnipsmit bas [sex efs snesexdex OF Bisse SIS oo bas rr ii 


; F > 2 ~~ 
neee ed neo tt (di + 8) sedmun xeigmos sdz io att Saogmeo 
. di a at Se ae 
enotivusattedye sidstivue 16546 3% 
sd + «6 (56 = 2) 
' 1 ~ . a = ho = 

) ( } r be en ee tee eres ae (2) | 

: S ite 

i 4 (6 * 2) 


: is 
mxoi sis ro satevns ns esd ci kttw 


-sop yd bavetdos ed m6> noteesxqxs Bint to no tt so£4 amie 


26 LSSqce Sivow mods 31 .8976N £br9609 xsfoq oF pnt, ek 


sa 
[~b) > +90) ute ~*s itd + 8) | x = @)% lina 
i 4 


uy: 


(2-) TABOMATOAA = (di + 6) p Yo tnamupts = 6 ox ft 


s - 
@ 


ages oes 


+2 


“ge + reN\ = |\{di + sii 


i 
... 


asenil betsesqez x10t s5ay OF exuheaaee xsiimie A 


he 


oijsibsup bestseqa: pailbasd ot ‘pewollot os sem 


teyisns efit at emeiderg siz to a . 


sneoes oite2zbssp 6 yd ‘betbasd ed 69 amaze a Loxtac 1D 
: a 

sebbnas od ” iw eas esis Laser? Rad bod segs 
sued ach eel hig Pye) Are ps Aili 
. - 4 : i i r. r t an) 7 
as ih ao ‘ - , 


Pr ee 


~ 


7 


49 
The original Laplace transform can be expanded as 


) 
Bian (s} a As + B a Cag top : 


Poesy se be ee ee ot tee eine 


op EES) (4-18) 


The constants A, Be C, and D can be solved for by multi- 
plying both sides of this equation by the square of the 
quadratic factor, differentiating with respect to s, and 
letting s + a + ib. The term F(t) corresponding to the 


repeated quadratic factor can be written in the form 


at 


F(t) = ae e®t tig(a + ib) | [SIN (bt + ©) = bt COS (bt + 0)] 


2b 
- bl ¢'(a + ib)| cos (bt + ©')|} + H(t) 


where ©' = arg [¢' (a + ib)] (4-19) 


3. Computer Implementation 


In the previous section the theory of the method 
of partial fractions for calculating the inverse Laplace 
transform of a function F(s) was presented. The applica- 
tion of this theory in developing a program to solve for 
the time response from a function polynomial presentation 
in the s domain will now be discussed. 

The transform function F(s) is a ratio of two poly- 
nomials P(s) and Q(s) as described in the previous section. 
In formulating a program to calculate the time response it 


was necessary to factor Q(s) into a set of linear and/or 
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quadratic terms. These terms or roots of Q(s) were cal- 
culated using the Newton-Bairstow method that is pro- 
grammed into subroutine SOLUN. 

The program that was developed to perform the in- 
verse Laplace transformation will handle up to twelve 
linear factors each appearing once, twice or three times 
and up to eight quadratic factors each appearing once or 
twice. This is an arbitrary limitation restricted only 
by the dimensions of the storage arrays. The program cal- 
culates the coefficients for the four general cases as 
given in Equations 4-7, 4-9, 4-13 and 4-18. This infor- 
mation is then used by another subroutine described in 
the next Section to calculate the time response. 

There are two vectors of data associated with the 
inverse transform, one for the linear roots and the other 
for the quadratic terms. For the linear terms it is 
necessary to calculate Aor A, and Aj or Ayr Ay and A, of 
Equation (4-9) depending on the number of times the linear 
root is repeated. The calculation of the time response 
contribution made by linear roots can be simplified by 
using a modified version of Equation 4-11. Instead the 
following expression for the constants, shown in Equation 


4-20, is used 


a: (a) (4-20) 


fu 


7h 


- - 
_ -{59.9xew (e)Q. to etoor to emiss SeedT ain was Seshiinia de 
-o1q eat ted? Horitem wosetisd&-noswe ad? partes betes a 


wuroe antodediae od mit melt 


-ti eft mxrotieq oF begolevsb esw 36d mexpoxrq SAT™ 
; ; oe 
sviews oF qu slbreér | lois amiotansyat 90 alqed 2219V 
aanist saudd tO ediws ,.s2ne izs9qqs dipss ee 
10 eDno pnixseqaqs doses etorosi oistsubsep sipie of 
rine betoiises toiste+imt textidys os at eid? 
-[60 mSxpoxrq siT .sYyBtts epstosye sds: 10 enotanemib 
an eseso fatsnep avet ss t etastottisoo sz e069 
-~tyoinik efit .8i-b bre ti- ,e-b Vb ey aan ad novi 
ni beditoesh enisuordse izedsons yd beeu asdt et aOLS 
~ .e: = 
-2enogeset emit sit ostsiluolso oF doLtgose Fxeq ong 
a eee = a 
siz dititw bestsiooegs sseb Ic NPS HS ows SiS SiSent ; 
. ae 
tsdjo adt bis vstoor teenmti sit uot sno  myotensxs seteval 
. “ A 
2i +i emtet xwsenil saz soO% .ensrxed oS $stbsup ent % o2 
r = 
to .A bas ,A , A 10 .,A bas (A ,.A stetuolso og yaseas 2 


1se9aiiL sas asmis to sedmun 
Senoges1 omit sd? to no.it 
yd Ssiitiqmfe ed mad esoor 
eit BestenI .{l+8 norisupa 
fioitsupa at 


‘ _ a. 
sit no pmibmegeb (&-b) aols supa 


aiwode ,etnsteno> sit x02 noleeexqne ba 


£ 


-bsssoqez ai 3ogx 
% 


L yd ebsam noistudin: 
‘to aotaiey beritbom s ening 


sivolso edt 


1S68nL! 


Sk 


where n = number of times linear root (s - a) is 
repeated. 
r = number of the coefficient to a maximum of n 


(Equation 4-9). 


The time response function that results when the AL. con- 


stants. are defined in this manner is shown below. 


at 


— 2 -_ 
E Cepatiehs care + MeN, thea ADE) (4-21) 


In determining the Ay and A, constants using 
Equation 4-20, it is necessary to calculate the first and 
second derivatives of ¢(s) and evaluate these derivatives 
at s =a. As stated before this function ¢(s) could be a 
ratio of two polynomials in a form such as 9(s) = N(s)/D(s). 
By using the relationship for the derivative of a function 


u(s)/v(s), simplified forms for the derivatives can be 


formed as shown below. 


sy) = Bis) NO) (s) - wis) pts) 
D(s) D(s) 


(2) (gy - (ls) - 2) D(s) N\ (s) = wis) p'?? (s) 
D(s) D(s) D(s 

The evaluation of the coefficients associated with 
a quadratic factor becomes slightly more difficult because 
the function must be evaluated for a complex number 


s=atiib. For a non-repeated quadratic factor there 
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are four parameters necessary to completely define the 


time response as can be seen from Equation 4-22. 


RE) oC e°" SIN tht + 6) (4-22) 


which is a simplified version of Equation 4-17 where 
C = |¢(a + ib)|/b. For the case of a quadratic factor 
that is repeated once, Equation 4-23 shows that six para- 


meters must be determined. This equation 


F(t) =ce®> {c' (SIN: (bt + 9) - bt COS (bt + 0)) 


tC. COS) tet. 104.) 3 (4-23) 


is the same as Equation 4-19 with some of the terms collect- 


ed and given the following meaning 


ed, : 


| ¢ (a + ib) |/2b 


o) 
Il 


-|¢'(a + ib) |/2b7 


4, Program Descriptions 


The first step in performing the inversion from 
the "s" domain to the time domain is the calculation of 
the roots of the denominator of the overall transfer func- 
tion of the system. A subroutine CALRT performs these 
calculations and begins by normalizing the coefficients of 
the denominator polynomial so the coefficient of the high- 
est power of s is unity. It then calculates the linear 


and complex roots of the normalized polynomial by using 
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the Newton-Bairstow technique which was programmed into 
subroutine SOLUN. 

After solving for the linear and quadratic roots, 

the coefficients of the inverse transform associated with 
each of these roots is calculated by subroutine INPFM. 
The coefficients for the linear roots are a, i oat A, and 
A. and for the non-repeated quadratic roots a, b, c and 
© while the repeated quadratic factor requires six para- 
meters a, b, -c, 9, C" and 6"; These coefficients were 
defined by Equations 4-20, 4-22, and 4-23. 

Since many digital computers including the IBM 1800 
do not have the complex arithmetic capability which is 
necessary for solving for the coefficients associated a 
complex root, three other subroutines were written for this 
function. They are PCMPY - to multiply two complex numbers 
together, PCDIV - to divide two complex numbers, and PCVAL - 


to evaluate a polynomial with real coefficients for a com- 


plex argument. 


(om Forcing Functions 


The main reason for calculating the closed loop trans- 
fer function is to allow the design engineer to calculate 
the effects that changes in one variable will have on the 
others in a unit or in a larger system made up of several 
component units. The changes of interest will generally be 
random variations - time functions of no special mathemati- 


cal shape - such as outdoor temperature variations, or feed 
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stream purity or flow fluctuations. 

Much valuable information about the effects of 
such changes can be obtained by considering the effects 
of a few standard changes. Some input functions which are 
frequently used for investigating the characteristics of 
a control system are the step, pulse, sinusoidal, ramp, or 
the impulse functions. 

The procedure is to pick some property of the sys- 
tem which is subject to change directly by forces outside 
the system, and treat this as an independent variable or 
input. Some other property of the system will be the de- 
pendent variable of interest and may be considered the 
output. Once a transfer function F(s) is calculated that 
shows the relationship between these two variables, a dis- 
turbance can be introduced at the input to see its effect 
on the output. This is said to be forcing the system and 
the five standard input functions mentioned above are called 
tne LOLrcing Lunctions, 

Three variables, namely, the setpoint, load one, or 
load two, can be used as the input or study variable while 
the output variable can be either C(t) or B(t) as defined 
in Figure 3-2. The subroutine COFF requests the user to 
enter the variable to be studied, the type of forcing 
function to be applied, and any parameters associated with 
the forcing function such as step size, or pulse height and 
width. The choice of the output variable was previously 


obtained by subroutine CLTYP as discussed in Section 6.3 
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of Chapter 3. Once all this data is successfully entered, 
the subroutine OCLTF calculates an overall transfer func- 
tion which may then be transferred to the programs used 
for calculating the inverse Laplace transformation and 


subsequently the time response. 


D. Time Response of the System 


1. Caiculation of the Transient Response Data 


The partial fraction coefficients must have been 
previously calculated before this Section of the program 
can proceed to calculate the time response data. An indi- 
cator is set any time the block diagram transfer functions 
are modified successfully. It is checked upon entering the 
time response calculation section of, the system and causes 
an error message being outputted to the user if it was not 
reset. It is reset after the partial fraction coefficients 
are calculated. 

The user is asked to supply data about the time 
response plot such as the end time, the number of data 
points to be calculated, and if it is an original plot or 
a replot of existing data. The time response will always 
be calculated from time t = 0 to the end time entered by 
the user if a new plot is requested. 

The time response data is calculated from the par- 
tial fraction coefficients by subroutine STIME. It does 
this by a summation of the contributions of each linear and 


quadratic factor at every time interval. Up to a maximum 
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of five hundred time intervals are allowed by the pro- 
gram. The response data is stored into disk file imme- 
diately after calculation. The replot option allows the 
user to display this time response data later without re- 
Calculating it. 

2. Display of the Transient Response Data 

After the time response data has been calculated 
and stored onto disk, the subroutine DSPLT requests the 
user to enter additional information about the display 
desired. This input includes which display device is to 
be used, scope or plotter, the plot direction, and the 
title of the graph. There are four possible graph types 
that can be obtained. They are as listed below and repre- 


sent the time response of the third order system 


G(s) H(s) = K /(s +.) (1/2 s+ UW) (1/73 s + 1) 


as Ka is varied and a unit step forcing function is applied 


to the setpoint. 


1) T}Vertical. plot (8% x ll inches), Figure 4-1 
ii) |}Horizontal plot (11 x 8% inches), Figure 4-2 
ii3) |} Top Half: plot (8% x 5% inches), Figure 4-3 


iv) Bottom Half plot (8% x 5% inches), Figure 4-3 


3. Performance Criteria 
The basic goal of a control system design is meeting 
performance specifications concerning the behavior of the 


system. Linear systems can be characterized by their 
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response to aperiodic input functions. Consequently, 
the design of linear systems is often specified in 
terms of the transient response for simple aperiodic 
inputs. Very often a step function input is used and 
some of the commonly used time domain specifications 
associated with the step function are presented below. 

3.1 SteadyState Calculations 

The calculation of the steady state value SSTAT 
of the transient response has been implemented so two 
alternative methods are available to the subroutine 
SSCAL for determining its value. The first and simplest 
approach calculates the slope of the last few time res- 
‘ponse data values. If this slope is small, steady state 
is assumed to have been attained and an average value 
of the last five data values is used for the steady 
state value. 

When the slope is not small, the second approach 
must be used. It results from the fact that complex 
conjugate roots of the characteristic equation give an 
exponentially damped sinusoidal transient response. An 
analysis showed that this was true for systems with a 
combination of real and complex roots. It was then sta- 
ted that the curves SSTAT tA eo would comprise an enve- 
lope which would always contain the time solution. The 
subroutine SSCAL uses a least squares method to determine 
these unknowns using the peak values of the responses. 


This envelope has been drawn onto Figure 4-4. 
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3.2 .Transient Response Criteria 


The transient response is normally analyzed by 
using a unit-step input as the forcing function. The 
response can then be characterized by the following: 

i) Percentage Overshoot - the percentage 
overshoot is found from the amount that the tran- 
sient response to a step input exceeds the final or 
steady state value of the output. It is usually 
expressed as a percentage of the final value and is 
a measure of relative stability. Allowable values 
for the peak overshoot range up to 25 or 30 percent 
if the settling time is fast. Overshoots higher 
than this are seldom acceptable. 

ii) Peak Time - the time required by the normal 
response or output to reach its first maximum. 

iii) Rise Time - the rise time is defined (26) as 
the time required for the response to a unit step 
function input to change from 10 to 90 percent of 
its final value. It gives a measure of the system 
response speed. 

Also included as one of the time domain spe- 
cifications and similar to the rise time is the time 
required for the normal response to the step input 
to reach its final value for the first time. 

iv) Delay Time - the usual definition of the time 
delay is that it is the time interval between the 


excitation and a measurable response. As a time 
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domain specification it is often defined as the 
time required for the response to a unit step to 
reach fifty percent of its final value and this 
is the time calculated by the program. 

v) Decay Ratio - this is the ratio of succes- 
Sive overshoots in the time response. 

vi) Settling Time - the time for the response to 
reach and remain within a specified percentage of 
the final value after a step forcing function is 
applied. Values of two to five per cent are often 
quoted as the specified percentages. The most 
commonly used value is five per cent so the program 
has been written using this value. 


Avene? COs Integral Criteria 


The purpose of a feedback control system is to mini- 


mize the output error as a function of time after some up- 


set has been introduced into the system. Both the magnitude 


of the error and the time over which an error exists should 


contribute to the definition of optimum control. Three 


of the most common error-integral criterions of performance 


are discussed below and are available on the system. 


i) IAE - Integral of the Absolute Value of the 
Error 


IAE = | Je(t)| dt 
1@) 


This criterion increases for either positive 


or negative errors. It is more sensitive to small errors 


than the ISE, 
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it) “ISH —“Integral "ef the “Squared Error 


ISE = | a ede 
O 


Both positive and negative values of the 
error increase the size of the integral. The 
ISE criterion usually yield responses that have 
small overshoots but long line-out times since 
small errors contribute little to the integral. 
iii) ITAE - Integral of Time Multiplied by 


Absolute Error. 
ITAE = | tle (et) |dt 
re) 


This criterion is insensitive to the initial 

errors but weights heavily errors occurring late in 

time. This means that the response will show a 

shorter total response time with less oscillation 

than the other criteria but it may exhibit larger 

overshoots initially. 

2.45 Display (ofithe Performance Criteria 

Because the performance criteria that were discussed 
above were so closely associated with the time response of 
the system, it was decided that it should be displayed to 
the user with the transient response. The programs that 
calculate these criteria allow the user to show the results 
on the scope or plotter as the bottom plot of a two figure 


plot as shown in Figure 4-4. Another option allows the 
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performance criteria to be typed out on the user's con- 
sole. This might be desirable when hard copy is required 


but plots of the time response have already been obtained 


or are not required. 
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CHAPTER V 


FREQUENCY RESPONSE 


A. Introduction 


Although the time response of a system to a func- 
tion of time impressed upon its input terminal is of 
basic importance much can be learned from the frequency 
domain analvsis. The essential feature of the frequency- 
domain method is that the description of the system is 
given in terms of its response to a diac. aldo varying 
input signal. 

The steady state response of a linear system to a 
sinusoidal input is also a sinusoidal of the same frequency. 
Therefore, the effect of the system upon the signal trans- 
mitted through it can be described in terms of two quanti- 
ties: 

a) The ratio of the amplitude of the response 
Sinusoid to the amplitude of the input sinu- 
soid, and 

b) The phase of the response sinusoid with res- 
pect to the input sinusoid. 

The amplitude ratio is termed ay gain of the sys- 
tem, and the relative phase angle, positive if the response 
leads the input and negative if the response lags the in- 
put, is the phase angle of the system. 

An important feature of the frequency response method 


is that the transfer function describing the sinusoidal 
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steady-state behavior of the system can be obtained from 
the transfer function simply by replacing the Laplace 
Operator s by Jw. The transfer function is thus a func- 
tion of complex variables, and, in general, can be repre- 
sented by the gain and the phase angle, each a function of 
the frequency. The complete frequency response is there- 
forerevealed by presenting both gain and phase angle over 
the frequency rangeO < w < ~. 

Frequency response methods are based on the open- 
loop transfer function G(s) H(s) as given in Figure 3-l. 
The plots are made by replacing s by jw in the transfer 
function, G(s) H(s). The performance of the closed-loop 
system is predicted from these plots. 

The following method of plotting the transfer func- 
tion G(jw) H(jw) are the most useful in the frequency res- 
ponse analysis and design of feedback control systems: 

a) Bode Plot (Corner Plot): A plot of the magni- 
tude (in decibels) versus the log of the 
frequency on rectangular (or semilog) coordi- 
nates and another plot of the phase angle 
versus the log of frequency also on either rec- 
tangular or semilog coordinates. 

b) Nyquist Plot (Polar Plot): A plot of the mag- 
nitude versus phase angle on polar coordinates 
as the frequency is varied from zero to infini- 
CY: 

c) Log-Modulus Plot (Gain Phase Plot): A plot of 


the magnitude (in decibels) versus the phase 
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angle of G(jw) H(jw) on rectangular coordi- 

nates with frequency as a varying parameter on 

the. curve. 
These graphs and the procedure by which they were obtained 
will be discussed in further detail in the following Sec- 


tions. 


B. The Bode Plot (Corner Plot) 


The frequency response of the system can be represent- 
ed by a plot of the magnitude (gain) versus the frequency 
and a plot of the phase angle versus the frequency. The 
Bode plot is also known sometimes as a corner plot or sim- 
ply the logarithmic plot of the open loop transfer function. 

The unit of gain most often used in Bode plots is 
the decibel which is defined as twenty times the common 
logarithm to the base ten of the magnitude. To facilitate 
the presentation of the frequency response characteristics 
over a large range of frequency, a logarithmic scale is 
used on the frequency axis. The phase angle is also plotted 
using semi-logarithmic coordinates. 

Since the decibel is a logarithmic measure of the 
gain, the over-all gain in decibels of two components cas- 
caded in series can be found by adding the gains in deci- 
bels of the individual components. The fact that the 
effect of cascading components can be determined by graphi- 


cal addition in the Bode diagram, rather than by graphical 
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multiplication as in the use of the Nyquist plot is a 
Significant advantage. 

To illustrate this method of exhibiting frequency- 
response characteristics, a Bode plot is given in Figure 


5-1 for the following transfer function: 


. ; oe LO..(s. > £0) 
Bac hate) = s (s + 5)(s # 2) 


It should be noted that in addition to the Bode plot in 
this figure, the gain and phase margins and crossover 

frequencies have been included to present as much info- 
mation about the system as possible. An optional list- 
ing of the frequency response data can be made on the 

line printer either when it is calculated or replotted. 
This will not be illustrated here but an example of the 


output may be found in the User's Manual (Appendix A). 


( Lhe Nyquist Dragram (Polar Plot) 

The polar plot of a transfer function G(s) H(s) 
is a plot of the magnitude versus the phase angle G(jw) 
H(jw) on polar coordinates as w is varied from zero to 
infinity. From a mathematical viewpoint, it is regarded 
as a mapping of the positive half of the imaginary axis 
of the plane of the complex variable s onto the plane of 
the function G(jw) H(jw). 

The Nyquist diagram is obtained by expressing the 


frequency function at a given frequency as a vector in 
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polar form: 


G(jw) H(jw) = Kw) / $(w) 


The locus of the tips of these vectors as w varies from 
0 to » is the Nyquist diagram. 

The subroutine NYPLT allows the user a choice of 
scales for the Nyquist diagram. It can be either automatic 
where the computer chooses the starting and ending coordi- 
nates for both axes by the subroutine COSGR or the user can 
specify these, thereby limiting the graph to a particular 
region. An illustration of these two methods, respectively, 


is given in Figure 5-2 and 5-3. 


Ds Log-Modulus Plot 


As mentioned in the introduction of this chapter, 
the log-modulus plot is a plot of the transfer function G(s) 
H(s) in decibels versus its phase angle in degrees with 
frequency as a parameter on the curve. The main advantage 
of using this set of coordinates is that such a plot can 
be superimposed on the Nichols chart to determine the rela- 
tive stability and the frequency response of the closed- 
loop system. When the gain constant K of the transfer 
function G(jw) H(jw) is varied, the plot is raised or 
lowered vertically according to the logarithmic scale. 

The magnitude versus phase plots are usually obtain- 
ed by first plotting the Bode plot and then using the same 


data to plot the log-modulus plot. The program follows a 
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similar approach by calculating the data for a Bode 
diagram and then plotting the data onto the log- 
modulus plot by subroutine LMPLT. In addition to the 
log-modulus graph, some frequency response specifica- 
tions are also listed at the bottom of the plot. An 


example of this plot is shown in Figure 5-4. 


E, Frequency Domain Specifications 


The performance of a control system can be charac- 
terized in the frequency domain by several different 
criteria. A description of some of the common criteria 
follows: 

1. Bandwidth 

The bandwidth is defined as the frequency at 
which the magnitude of the closed loop transfer function 
M(jw) has dropped to 70.7 per cent of its zero frequency 
level, or 3 db. down from the zero frequency level. 

The bandwidth gives the range of frequencies for 
which the system gain is adequate to insure good trans- 
mission of the signal. A large bandwidth usually indi- 
cates that higher frequency signals will be passed on to 
the output. Thus, the transient response may have a 
faster rise time accompanied by a larger overshoot. Con- 
versely, if the band width is small, only low-frequency 
Signals are passed; therefore, the time response will 


generally be slow and sluggish. 
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2. Peak Resonance, M 

This is defined as the maximum value of the closed 
loop transfer function M(jw) and also gives an indication 
of the relative stability of the system. The peak resonance 
has been shown to be related to peak overshoot for second 
order transfer functions (13, 22, 26) and any system deter- 
mined by a single pair of conjugate poles could be satis- 
factorily designed using the M, as a specification. It 
is apparent from Figure ee, given by Murphy (26) that a 
high M,, corresponds to a large overshoot in the time res- 
ponse. For design purposes, it is generally accepted that 
an optimum value of Se Lee ieietorl. > which results in an 
overshoot of from 15 to 30 per cent. 

3. Resonant Frequency, Wy 

This is the frequency at which the peak resonance 
M, occurs. 

4. Gain Margin 

The gain margin is a measure of the closeness of 
the phase-crossover point to the critical point. The 
gain margin frequency is first calculated as the frequency 
at which the phase crossover occurs (phase angle equals - 
180 degrees) and then the gain margin is found by deter- 
mining the magnitude in decibels at the gain margin fre- 
quency. 

The physical significance of gain margin is usually 
that it is the amount of gain in decibels that can be 


allowed to increase in the loop before the system reaches 


ay PP i 
~~ 4 : 
e ? “ie 
i on 
“oS ano As gq Sin 
_M , S2oat poEer 2 ‘= 
5 | *. - J 
a ] 
. = P ” ; 
beszolo 3 ‘Oo eylsay munmixem as eS jualin ati atnT ¥ 
; an 7 as - 
ne t ms a@evip offs bas (.f)M aossonpe ustaasxt gool 
é a e 9 
a lies eel 
IDI q ort teaseve to yotiftdete eviteles es 
bacoose z0ot 20 “evo ABssq Of bosslet ed od mworie: need | 
i = ill — 5. . c helat a fal > > 3 
2 G6 quar Dns - ‘= 5 ~ i} wo LaIsnwvt 1928ns7 
Is i bfuonp sxefoq stsputnom to uisq eipntia s yd 
1rOLsAD oeqe 6 8&5 _M afa potas Beapiseb ylt 
sari? (88) vaqieM yd nev .S-@ suveid mould ee 
win tm 
“231% Sms Snivy iu rOeHhSIsSsvo spas B&B O5 BDNOGesIi0oD Pf. > = 
tedt botqesos yilsxensp = t ,eseoqxdg aptesb 10% ~,: 
Is ai esives foftriw-d@.f og 1.f£ @£ a” to apisv 
~ wid 4 
tie> seq Of ot Gf Moz? to 
4" Xpesuy pHi tmemnozet «f° 
 < A TT TS IY 
soasnrotet Assgq sis doidw ts yousnpex3 sis ai eiAT 
| , 2% 
aipasM atsD .t 
20 eBesnseo sit to estuesenm s ef sipasm aiep onT, 
‘. 
eff .tnteq isorsixuo eda oF taieg- revoaeez-sasdq ‘od 3 
ms 
yomeppet? sds 85 Hhessiuolso sexirri ei Ys apopert okgaa Fl 
~ afseups sipn sig) exo: “WsvoeBZorD: seetg ois doidw 9: 


—~teteb vd bavot 


-S1i nipism niep aft 36 


we a 


> 081 


bins (2oanpe ey 


77 


instability. The gain margin must be positive when ex- 
pressed in decibels for a stable system. 

5. Phase Margin 

The phase margin is a measure of the closeness 
of. the gain-crossover point to the critical point. It 
is calculated by finding the phase angle of G(Jjw) H(jw) 
at the gain crossover frequency (point where the magni- 
tude is 0 decibels) and then adding 180 degrees. The 
phase margin is the amount of phase shift at the gain 
crossover frequency that would just produce instability. 
The phase margin must be positive for a stable system 
and a commonly used value for design purposes is 45 
degrees. 

It is usually preferable to evaluate the gain and 
phase margins of a control system from its Bode plot. 
The reason is simply that the Bode plot is easy to con- 
struct, and the margins can be obtained by inspection. 
Because this was true for the cases when the Bode plot 
was drawn by hand, the frequency analysis programs cal- 
culate the gain margin and phase margin and display the 
values on the Bode plot. As an illustrative example, 
consider the open-loop transfer function of a unity feed- 


back system. 
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The Bode plot is shown in Figure 5-5 with the log- 
modulus plot (Figure 5-6) giving the additional speci- 


fications of peak resonance and resonant frequency. 


78 


- ¥ a ‘ i 
% ; o 
i - ips 
ol ods dplw 2- i mc mwosea at tolq sbo | 
: ca - 
: 7 
- notszibbs ene ¢ o (8-2@-saupiS) solgq evlubom 


nenoaes Jmeq to anciteott 


20 SL ll Dn rl Di 


MAGNITUDE 


-149» 9000 


C 


Bebe pet ed Ss 


ANGLE 


PHASE 
-275+0000 


20» 9000 


mo OOO 


79 


BOLE DLAGRAM 


ee ee eee tt A) A |) 


GAIN MARGIN = .-34-807 DECIBELS 
GAIN CROSSOVER =) 0-S81b CO RAD-/SEL- * 


Sr" Sus A oan ESE Da eR GL LAA ceca ee a Ma ed SE aoa tH 


PHASE MARGIN = 77+774 DEGREES cae 
PHASE CRESSLVER © +158 C2 RAD-/SEL- a ae 


Sig ce) a> Gal ene Lean meena et ae jeepers 
6*41000 S89 «Soe 
bg Aa gee eA Le es eer 


FIGURE 5-5 Second Bode Diagram 


~~ ¢ 7 __ 
i ‘ 6 | . _ 
; ' 4Puhta) (ett Smee! | 
an é . } ) o ; ¥ ts 
| | nie |. 
re - — 7 
os a 


= 


= 


be. 1 rar : 


= 


1] oa 


Wiehe eee BE ees 


LOG; - MODULUS OTAGRAM 


80 


KUO 2€ REF: 22 ; PAGES 401-404 ) 

3 } 

. F 

A 2 

4h a 
peg! 
eee 
sai 
thie 
va 
a 
r 
vo 
# 
af 
af 
£ 

O 

i) 

a 

3 

= + feo eel es a 

-275 +0008 =§5«°GO0CC 

a eee aie =n RES 

PERFORMANCE CRITERIA FOR ABGVE DSTAGRAM 
Prinoe URJSSEVER = oe eel eo Se 
~ (SAIN MARGIN ehhices AnijOdoublich [BELLS 
GAIN CROSSOVER = Cia te AL ae aie 
PHASE MARGIN = PAK Yiog AAEM IS ES 
PEAK RESONANCE me eee Cee 
RESONANT FREQUENCY = Ce] GOOG PAL ae -“Seiee 


FIGURE 5-6 Second Log-Modulus Diagram 


81 
CHAPTER VI 


ROOT LOCUS ANALYSIS 


A. Introduction 


The root locus diagram was introduced by W.R. 
Evans (14) in 1948. This technique has been improved 
and applied to the analysis and design of control sys- 
tems in recent years. 

The root locus is defined to be the trajectory 
in the s-plane followed by the poles of the closed-loop 
transfer function as some parameter of the corresponding 
system is varied continuously. The characteristic equa- 


tion of any system can be written in the form 


1 + KG(s) H(s) = 0 (6-1) 


The function G(s) H(s) in general may be a rational poly- 
nomial and/or a transcendental function of s. The graphi- 


cal methods of Evans based on the angle condition 


/Gisoiy Hitsdp f= 4744) 2ntie nyse 04 ty) 202 Dia ee 4652) 


has been well developed, but it is a trial and error 
procedure. It can therefore be tedious and inaccurate 
and in such cases, an analytical solution is preferable. 


B. Mathematical Basis of Calculation Methods 


The open loop transfer function G(s) H(s) in 
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Equation 6-1 can be written as a ratio of polynomials 
N(s) and D(s) of degrees m and n, respectively, where 
m < n and K is a parameter which is proportional to the 


system gain 


Gis) H(s) = HS 
m m-1 m-2 
P Seitigaenj® +a.» s tooo. + ay 
n my en 2 
Sg aes mpi se *% b-2 Ss + 5 a eh 
(6-3) 


The complete frequency s + jw can now be substitu- 
ted for the variable s, and by separating the real and 
imaginary parts of G(c + jw) H(o + jw), the equations 


governing the root locus are given by 


Hy] 
1 


Re G(o + jw) H(o + jw) (6-4) 


Oo AIR 


Im G(o + jw) H(o + jw) (6-5) 


It was stated in Equation 6-3 that G(s) H(s) could be given 
as the ratio of two polynomials N(s)/D(s), this substitu- 


tion can now be made into Equations 6-4 and 6-5. 


Re N(o + jw) Re D(o + jw) + Im N(o + jw) Im D(o + Jw) 
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Re Dic + jw) Im N(o + jw) —- Re N(o + ju) Im D(o + jw) = 0 
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Equation 6-6 is called the "gain equation" and Equation 
6-7 is referred to as the "root locus equation". Ifa 


pole of the open loop transfer function, o, + jw, satis- 


1 
fies Equation 6-7, then the value of the gain correspond- 
ing to this pole can be obtained from Equation 6-6. 

If N(s) and D(s) are analytic, say a polynomial 


in s (with real coefficients), they can be expanded into 


a Taylor series about any point (0, jw): 


. 2 3 
N (o 55 jw) == Nic) “+ sun 2) (co) a ow?) (a) is 24 n (3) (co) a eA 
| (1) THe AES Tere. 
D(o + 47m) = Wo )ichyjwd (s) - 5T D {o) - i D Ca ite... 
(6-8) 
where y 4) = the ith derivative with respect to o. 


- Equation 6-8 can be substituted into the Equations 6-6 and 
6-7 and after some simplifications new expressions can be 
obtained for describing the root locus. 


The root locus equation can be expressed as 
[R,(c) - w°R,(c) + w Re(c) # ..-] = 0 (6-9) 
ae 3 5 - eoe 


while the gain equation is given by 


Roo) - w?R, (0) + wR, (o) - wR, (a) 
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The coefficients of the powers of ae in the previous two 
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k (xr) (ker) 
#8 “7,2 N (a) “5 (0) a 
Ry (9) s iis (~1) a k-r)! (6-11) 
k (r) (k-r) 
is ee ane (sc) D (oc) = 
Stak) 4. 45us _ oon (6-12) 


The degrees of Equation 6-11 and 6-12 are given by 


(n + m- k) and (2n - k), respectively. 


A number of observations can be made concerning the 


root locus equation. 


a) 


b) 


Cc) 


d) 


e) 


f) 


The real axis w = 0 is always on the root locus. 
The locus off the axis is expressed by a finite 
series in A 

The last term in Equation 6-9 is given by 


(-z) (ntm-1)/2  ntm-1 Ragm(0) if ntm is odd, and 


ey 2 ate Rigm<1 6?) tf nt+m is even. 


The root locus equation represents an algebraic 
curve and is derived for a given ratio of poly- 
nomials without any prior knowledge of the loca- 
tion of the poles and zeros of the open loop 
transfer function. 

Terms of powers higher than fs are present in 
Equation (6-9) only if the sum of the degrees of 
the numerator and denominator of the open loop 
transfer function exceeeds 6. 

The breakaway points on the real axis are the 


intersections of the root locus and the real 
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axis and are given by the solutions to the 


equation 


R, (9) = 0 


The coefficients of the powers of ae in the Taylor 
series expansion of N(s) and D(s) given by Equation 6-8 


are first arranged in an array called the root locus array. 


Ry (0) R, (9) Ro (9) R3 (>) Ri 6%) R69) 
(2) C3) eas eee) 
N16) | WS!) pepl: NoculP gn iealale phen( Mi gEza(eh) aub- 0 
; (2) (3) (m) (n) 
D(o) pg) D — D = coe D ee eee D wl 


‘In this array the odd order differentiations of the numerator 
N(t) are negative. The coefficients of R, (o) are obtained by 
taking the sum of the (k + 1) cross products. For example 


the calculation of Ro (9) would be as follows: 


Ry (co) = N(o) =n‘) (5) 2 a7 = 
ome 
D(c) D! (a) P = ) 
(2) (2) 
= N(o) PEAS) - w) (co) DM (o) + BS dic) 


The coefficients Dylods Do(o), D,(o), etc. are obtained in 


an analogous manner. 
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C. Computer Implementation of Root Locus Algorithm 


The root locus algorithm that is described in the 
previous section and summarized in the form of Equation 
(6-9) and (6-10) was implemented in a FORTRAN program. 

The derivatives m andn derivatives of the numerator and 
denominator polynomials, respectively, are first calcula- 
ted. This enables the root locus array to be calculated. 
Subroutine RCOEF determines the R, (9) coefficients by 
Equation (6-11) while subroutine DCOEF and Equation (6-12) 
are used to calculate the Dy (9) coefficients. 

The root locus equation is then evaluated by sub- 
routine RLEVL to determine all the values of the frequency 
that satisfy the equation. The real axis is divided into 
fifty segments and its length is governed by the low and 
aah frequency limits that the user was requested by sub- 
routine RLINP to enter. Successively these fifty values 
are substituted into the root locus equation and the solu- 
tions that give a positive value of the imaginary fre- 
quency component are determined. The frequency with a 
zero imaginary component is always a solution. 

The values of the gain associated with each of the 
frequencies determined by RLEVL is now calculated from the 
gain Equation (6-10) by the subroutine GAIN. The root 
locus algorithm that has been used results in frequencies 
that give both positive and negative values of the gain 
parameter K which implies negative and positive feedback 


respectively. 
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87 
D. Display of the Root Locus Diagram 


The root locus data is displayed onto a point- 

plot by subroutine RLPLT. Both the positive and negative 
values of the gain are placed onto this plot. They are 
differentiated by plotting the frequencies points with 
positive gain as small plus signs (+) while those points 
with negative gain are plotted as small crosses (x). The 
gain values are not plotted onto the graph but can be ob- 
tained from the optional printout of the values used to 
plot theroot locus. Three examples of the root locus dia- 


grams obtained are given here as 


i) G(s) H(s) SEES) Figure 6-1 


s + 16s” + 108s + 400s + 800 
Figure 6-2 
ii?) G(s). H(s) Figure 6-3 
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FIGURE 6-1 Parst Example of Root Locus. Diagram 
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The preceding chapters have presented the mathemati- 


cal basis for the different computations, and given example 


plots and a brief description of each section of CSDAP. 


This chapter will illustrate one typical use of CSDAP, namely 


a study of the effect of changing controller types and con- 


stants on the output of a control system. The various plots 


are given in this chapter, while computer generated documen- 


tation is given in Appendix B. Detailed operating procedures 


and further examples of all possible information that can be 


obtained from the CSDAP system is given in Appendix A which 


has been written to serve as a User's Manual. 

The control system studied is a single 
loop consisting of a process described by two 
transfer functions with time constants of 1.0 
a feedback element described as a first order 
L1OnW With a time constant of 0.33. The block 


illustrated below. 
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This control loop will be used to illustrate how 
the CSDAP system can be used to study the effect of 
different controllers on the response of the system C(t) 
when a unit step input is applied to the system as a load 
disturbance, (L1(t)). The results will be expressed by 
means of the transient response, the frequency response 
(Bode, Nyquist and Log Modulus), and the root locus dia- 
grams. 

The CSDAP system is initiated as a non-process 
program from the card reader. A listing of the control 
digit and data switch functions is obtained on the type- 
writer by previously setting data switch zero in the "ON" 
position (see section A.2.2. of Appendix A). The first 
function that must be performed is the configuration of 
the control system. This is accomplished by entering a 
CONTROL DIGIT of 1 and then the correct Mdheeey functions. 
Data switch 2 is set to the "ON" position before entering 
the CONTROL DIGIT of 1 and causes a listing of the trans- 
fer functions to be made after they are correctly entered 
into the system. Since, for the initial configuration, a 
proportional controller with a gain of 1 is desired, it is 
not necessary to enter the control one transfer function 
because any transfer function that is not entered by the 
user is assumed to be equal to unity. After the process 
and feedback element transfer functions have both been 


entered the user terminates this section of the program by 
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use of the END-OF-FIELD (EOF) character. This results in 

the statement "ILLEGAL TRANSFER FUNCTION DATA" being typed 
out to the user and causes the program to request the next 
CONTROL DIGIT. 

The user now instructs the CSDAP system to calculate 
the closed loop transfer function and the inverse Laplace 
transformation for the system response desired. Data 
Switch 2 is reset to the "OFF" position and data switch 3 
is set to the "ON" position to give the input/output variable 
and forcing function selection options. A CONTROL DIGIT of 
2 is then entered. The system response C(t) to a unit step 
in the load variable (L1) is desired so the set of variables 
shown on page B-2 is entered. Data switch 3 is then reset 
to the "OFF" position. 

Since the user has a CSDAP User's Manual (Appendix A) 
to consult for variable options he resets all the data 
switches to the "OFF" position to restrict further assistance 
from the program. The user then proceeds to have CSDAP cal- 
culate the transient response of the configurated control 
system by entering 3 for the requested CONTROL DIGIT. A 
zero is entered for the time response replot digit since 
the transient response data has not been previously calcula- 
ted. On the basis of information provided by the user, 250 
transient response values over a time range from 0-25 seconds 
are calculated. The data is then plotted onto the storage 


oscilloscope as an upper-half plot with no title. The 
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performance criteria calculated from these transient res- 
ponse values are then displayed on the bottom half of this 
oscilloscope by entering a CONTROL DIGIT of 4 and an output 
device option of 1. After reviewing the oscilloscope dis- 
play a plot of the transient response and performance 
criteria was obtained on the digital plotter by specifying 
a plot device of 0 instead of 1 in the previous parameter 
specifications that caused the display to be put on the 
oscilloscope. 

The next step in the study was to vary the propor- 
tional gain constant to satisfy the arbitrarily chosen 
criteria of about 30 percent overshoot. The proportional 
gain was changed by entering a new transfer function for 
controller one, Cl(s). The gain was set at 0.5 and the 
transient response calculated and displayed on the oscillo- 
scope and the performance criteria printed on the type- 
writer. The copy of.the transient response curve that was 
plotted from the data calculated for the oscilloscope dis- 
play is shown in the upper-half of Figure 7-2. The control- 
ler gain was then set to 0.75 and the transient response 
curve that resulted is shown as the bottom-half of Figure 
7-2. The performance criteria were then displayed on the 
oscilloscope. It was necessary for the user to manually 
erase the oscilloscope before requesting the performance 
criteria to be displayed because normally the criteria are 
displayed below an upper-half page plot and the oscilloscope 


would have been erased for that plot. 
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TRANSIENT RESPONSE 
1-4, PROPORTIONAL CONTROL ONLY 


OG 5 10 ite) 20 rere 
TIME 
peep vane an Pherae ABOVE RESPONSE 
SEAL aS Wale = C-7500 
PERCENT OVERSHOOT = 47°18 
PEAK TIME - 41-2048 
es NE. sO. Le) toe) O-6370 
RISE TIME 10 TO SO = ievcace 
DELAY TIME = O- 3420 
BEE AY RATA = G2 iid 
4-51.60 


Git IL LINis* iM = 
BRRUR INUEGRAGS ~ 


TAL = OC eS ets me 
{SES O- 246/75 
Lee gal Veet ame oo 


Figure 7-1. Unity Gain Proportional Control 


OUTPUT VARTABLE 


See vAR LARUE 


96 
Pee NS LENT ee ee iNSE 


ea (ind | ee at = Oe SO) 


+ =} 
@) 5 10 oh) 2G C5 
TIME 
Ee ANS LEI ee ar GINSE: : 
4 1 cel ty 2 12 ed Dera as a gta 
8 
S 
2 
$$ +— ——+ 
O = 10 Abe eo oma 


TIM 
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The displayed performance criteria indicated an 
overshoot of about 35 percent. On the basis of the 
results a value of the gain constant equal to 0.65 was 
chosen and a plot of the transient response and the list 
of associated performance criteria obtained on the digital 
plotter. Before entering the CONTROL DIGIT 3, data switch 
6 was set to the "ON" position to obtain the tabulation 
of the transient response data on the line printer. These 
values are given on page B-22 of Appendix B and are the same 
as those used for the plot shown in Figure 7-3. After the 
listing had been completed the data switch was reset to the 
"OFF" position. 

The user then decided to list the overall transfer 
function for the existing system. This was done by enter- 
ing a CONTROL DIGiT of 7 and a list device option of 0 
indicating the typewriter. A transfer function name of 
TM was entered to obtain the overall transfer function 
listing (page A-23 ). The resulting transfer function list- 
ing is shown on page B-7. 

The controller gain constant was then changed back 
to unity and the frequency response and root locus diagrams 
obtained. The use of logarithmic scales in addition to the 
unity gain offers substantial advantages when making cal- 
culations about the system. The gain margin can be obtained 


direct from the Bode diagram since it is the difference 
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between the overall gain and unity at a frequency equal 
to 180 degrees lag. 

The frequency response data was calculated for a 
frequency range from 0.01 to 100.0 radians per second. 
After displaying the Bode diagram on the oscilloscope it 
was decided to set the range from 0.1 to 100.0 radians 
per second. This removed the first decade from the plot 
which was of a constant phase angle. The Nyquist and Log 
Modulus plots are constructed using the data previously 
calculated for the Bode diagram. An automatically scaled 
Nyquist diagram is first displayed on the oscilloscope 
before setting the range of vaiues for each axis of the 
plot to be drawn on the digital plotter. The Bode, 
Nyquist and Log Modulus diagrams are shown in Figures 7-4, 
7-5, and 7-6 respectively. After the desired frequency 
range is chosen for the root locus diagram by a similar 
method to that used for the Bode diagram it was plotted 
and shown as Figure 7-7. 

The next phase of the study was to investigate the 
effect on the system of a proportional - integral control- 
ler when the integral time constant was varied. The trans- 
ient response curves for integral times of 0.5 and 0.75 
and a Hroportional constant of 0.65 are given in Figure 
7-8. The transient response curve and the performance 
criteria for an integral time constant of 1.10 are shown 


in Figure 7-9. The value of the integral time equal to 


ar ‘ ’ ’ 
: i ' 2 - : 
; He ‘ce 
g - 
} | ye 
, BUpS yonsupestt £6 JB yIinu ONS ALep Lisisve edt vad NIS a 


' 
pad £ .aestpeb et. 
“s 
103 SetseivolsS asw s95b senoqees ywaorms x2 sift. a 


O3 
“ae 


18q ensibsy 0 ot 10.0 raat spist yous ps’ 
. P - ‘ * = a 
i sqove0llioeso sd3z no msypsib shoe sna oniyelgelb x04 


od 1.0 moti spensa oda Jos oF bebiseb) as 
a 
q it mowi sbheosb Jers 2i3 bevoms: efnT .Baopes teg 
a 


— 


pol Snes selupy eff .sipas eesdq snatemoo 6 to enw dosti 


te 
a5 W 


: ' 
; > ae - = z 

faveivetg sisb ons enste storswtanos s16 aszolg aulubow 
‘pelsos ylisor A perch ebod siz xzo2 betsivolss 


% 


; _. toe 
2>f[lLineo sit no SeysiqeLh texzrs es mexpsib _— upyl 


! % : a wk = Lean 
e. “ F, gc my uA { Hine < ele 


—— 


* : 4 4 —- an 

,f uspld af dwode ©16 smsipelS eolonpom pod Das teic : 
: 
amsipesai Ositasd eas AassiséA YA ovisoageox a-~ bas , c= t 
uslimie s yd msupbBib avool tooxr sad. 20% nesodo af Spmst 
73er 
betjclq asw ti maypralb sbof sdy yo? Besse Bans oF aoe te 2m 

i f - 


.T-) siupit es avo | Dns 
sit sispiszasval of asyw ybuse sas Fo szariq sxen oat 
_=Loutnoo ibtpsdni - [snolsicgorg 5 to nosaye <orl Fet-) ts 


“ener? oT Heltavy esw tnbstends suit L832 podad ont | hadw 1s 


potit. tok ai sanogee 


c 3 e o a ae 7 ba A 


¢ 
eonamiotisq eds bas puse~ bch : 


ik 


2.925 or. - M to” Insten >: ea ‘5 


Co 


ms 


giclee 


— 


MAGNITUDE 


-105 » GGOO 


eae 


4 


Sess 


Pak se  ANGEE 


~280 » O000 


415» Q000 


40 » 9000 


100 


BOLE LITAGRAM 
PROPORTIONAL CONTROL 


A 


GAIN MARGIN = -10-457 DECIBELS \ 
GAIN CROSSOVER = 0-170 O41 RAD-/SEL: 
tp tt 


PHASE MARGIN = S0°i106 CEG EES me 
ASE: CROSEDVER = =e He aote Ol. RAD SEG. ee ae 


Fora ee eae lode gh, 2s hed aR nail OTe AO Ra Pel sg 


O* 1000 99.9999 
Pi Oe ie ALLL AS Ae aN 


Figure 7-4. Bode Diagram for Proportional Control Only. 


IMAGINARY 


D» BOO0E O00 


ri ik | AY ciate OH 


BE 


NYQUIST DIAGRAM 
PROPORTIONAL CONTROL -- 


: 


ek i ; 
Sere rOOGH O41. Cs BOOS CL. 


ale 


Proures/—5. — Nyquist Diagram For Proportional Control 
Only. 


Mt 


= 


=F eivol 
=) 
: 


valde Mees 


e 


Pe 


~ 


E 


MAGNITUO 


LOG - MODULUS DIAGRAM 


PROPGRTIONAL CONTROL 


415» 0000 


- 105 » 0006 


-269 - COCO 
PHASE ANGLE 


- DEGREES 


PERFORMANCE CRITERTA FOR ABOVE 


Beaeee OROUSSCVER 
GAIN MARGIN 

GAIN LROSSOVER 
PHASE MARGIN 

PEAK RESORNANEE 
RESONANT FREQUENCY 


Peer kis tt hg 


=e 
i 
30 
= 


ee 
o 

Q 

g 
bing 
ao 


i ar 
Aaa) Se) 
PAG es 
1Ob@ 
cd oe) 


Ole ter 


102 


a pp 


AQ+00@ 


DLAGRAM 


Pere Sa (oe 
i ee 
Ale a 7 ile 
oe hs tees 

fad ke) ne 
RAL ® “GEC « 


Figure 7-6. lLog-Modulus Diagram for Porportional Control 


Only. 


ale 


one | De eh ne 


ti 


1 NE RE ee re Ps ess 


14° Q0G0 


at oa ©) 0) 018) 


103 


PLEA SCC IEING EREATERAM 
PROPORT IGNAL: CONTROL 


Poe roll Ve ore 
a NEGATIVE GAIN 


-8B *QO00 


Prguce y=. 


ee 


Pave X 1 


Root Locus Diagram For Proportional 
Cconcrol Only. 


fc 


> = 
a 


> 
Y 
) : 
") a AZ : 


- 
> 


} 
* 


‘ 
x “000 
a 


7 
eee 


——t ae Ne 
> 


fie 


z 


ie 


1 


ares 


» 


= a —_ a = ed 
- — <a a) " 
Saha a ee Re 


be VA LAS Ae 


nS 
= 
- 
roa 
Ss 


Goria VARIABLE 


ee er PEN INGE, 
Seer AN aie bee oles el Oe 


Peele Ne Reo INS 
eee aN Weis KE =Osb530lT = O45 


3 
5 
vi 
23 
my + 4 —+- + a 
0 5 10 15 20 25 
TIME. 


Figure 7-8. Variable Integral Action For Proportional 


iicedtadecOontro. . 


ae + 


: | ere Sn eS 


GISTs) OVA Acie 


105 


Pie ee Set NE 
een IN Re ets aC 


as 
“G 
3 
O 
ved tt rt 
Hg = 10 15 mS, Coy 
TIME. 
PERFORMANCE CRITERIA FOR ABOVE RESPONSE 
ariAlgyve SEALE: = O-Q000 
PERCENT OVERSHOOT = 0-00 
PEAK TIME = cez= 


ples he aA ee Oa LG al S Be 
Pub te Mi Mert) 9° 20 += 
EISAY wt IME. | aig. 
WERAY RATT S a 


ees ieee Co 
O 
‘eS 
O 
[Ss 


Sey iN -ToL Me. dete erat 
ERROR INTEGRALS 

TAGs uss 2-aehs4 

eRe EA on PAM tee 

ITAE = peg ose 4 


Figure 7-9. Proportional-Integral Control 


106 


1.10 was arbitrarily chosen as the best integral time 
constant. 

The transfer functions were then listed. A partial 
listing showing only the control one, process two, and 
feedback one transfer functions was made on the typewriter, 
and it is shown on page B-13. A complete listing of all 
the transfer functions was obtained on the line printer as 
shown on pages B-23 to B-25 of Appendix B. The Bode, 
Nyquist, Log Modulus and Root Locus diagrams obtained from 
the plotter by a-similar procedure to that used for the 
proportional control only are shown as Figures 7-10, 7-11, 
7-12 and 7-13, respectively. Before the final Bode dia- 
gram was plotted data switch 10 was set in the "ON" posi- 
tion to provide a listing of the frequency response values. 
The line printer listing is given on pages B-26 and B-27 
in Appendix B. The root locus values also listed to the 
line printer are given as pages B-28 and B-29 in Appendix 
B by having data switch 14 set in the "ON" position when the 
final root locus diagram was being plotted. 

The effect of adding derivative action to the pre- 
viously determined proportional integral controller was 
then studied. Figure 7-14 shows the transient response 
curves for values of the derivative constant equal to 0.5 
and 0.25 while Figure 7-15 shows the transient response 
curve and the performance criteria for a third value of the 


derivative constant that was employed. For this value of 
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the derivative constant equal to 0.01 the plots of the 
Bode diagram, Pee Ryantst diagram, the Log Modulus diagram 
and the root locus diagram are shown in Figure 7-16, 7-17, 
7-18, and 7-19 respectively. 

The program was then terminated by entering a CONTROL 


DIGIT of 9. 
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Figure 7-16. Bode Diagram For Proportional-Integral- 
Derivative Control. 
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CHAPTER Vitt 


CONCLUSIONS 


The CONTROL SYSTEM DESIGN AND ANALYSIS. PROGRAM 
(CSDAP) presented allows the user to design and analyze 
continuous linear feedback control systems that can be 
represented in the familiar block diagram and transfer 
function form. This system provides the control engineer 
with an easily used program that operates in a time shar- 
ing mode on a medium sized computer. Most conventional 
control systems design and analysis techniques using the 
transient response, the frequency response, and the root 
locus diagrams have been implemented in CSDAP. 

The examples given in Chapter VII, Appendices A and 
B, and throughout some of the other chapters confirm that 
the original system design specifications for the CSDAP 
system have been accomplished. Some of the specific features 
of CSDAP that should be noted are the following: 

(i) The program provides for completely documented 
problems in the form of a typewriter listing of 
the user program correspondence, calculated 
values available in the form of tables on the 
line printer, and graphs of the response curves 
on the Calcomp digital plotter. 

(ii) The program itself provides the user with a con- 


Siderable amount of assistance as data switch 
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(iv) 
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options in all sections of the program. The user 
can gradually reduce this assistance as he acquires 
more experience with the use of the CSDAP system. 
Individual block transfer functions can be modified 
without re-entering the whole control system con- 
figuration. 

The user is provided with the ability to list all 
Or individual transfer functions to the typewriter 
or line printer as he desires. 

The fact that the program has been developed in 
sections allows the user to alter problem specifi- 
cations and parameters before rerunning the system 
analysis. For example, this means that the user 
can respecify the desired axis co-ordinates for a 
Nyquist diagram without recalculating the frequency 
response data values. 

The effect on the control system behavior of a 
number of forcing functions when applied to the 
input variable of the system can be studied. 
Forcing functions that are presently implemented 
include the step, pulse, sine wave, ramp and 
impulse. Possible input variables to the system 
include the setpoint and the two load variables. 
Transient response values can be calculated from 


the open loop transfer function or from the closed 
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loop transfer function when the output 

variable is taken as a signal either before or 
after the feedback element of the outer control 
loop. 

The program allows a time delay to be included 
in either the forward or the feedback path of 
the control loop by using the Pade and Taylor 
series expansion approximation methods. The 
transient response data can be displayed in the 
form of a table of values on the line printer 
or as a plot on either the digital plotter or 
the storage oscilloscope. 

A number of performance criteria such as the 
steady-state value, percent overshoot, peak 
time, rise time, delay time, decay ratio, 
settling time, and the error integrals associa- 
ted with the transient response data can be 
calcuiated and displayed on the digital plot- 
ter, the storage oscilloscope or the typewriter. 
The program generates the frequency response 
data and displays it in the form of a Bode, 
Nyquist, Log Modulus or Root Locus diagram or 
lists it in the form of a table of values on 


the line printer. 
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(xi) The frequency response criteria can be calcu- 

lated from the data stored on the disk. The 
gain and phase margins and crossover values are 
displayed on the Bode diagram. In addition to 
these values, the peak resonance and the reson- 
ant frequency are listed on the Log Modulus plot. 

(xii) Options are present for replotting previously 
calculated transient response, frequency response, 
and root locus data that was stored in the disk 
files. The structure of these data files is 
described in the Systems Manual and should allow 
other programs to calculate and store the data in 


these files for use by the CSDAP system. 
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CHAPTER IX 


RECOMMENDATIONS 


Experience gained while using the program to gener- 
ate the examples shown in this thesis has suggested a 
number of changes and additions that might be desirable 
for better program performance. 

The first modification is concerned with the entry 
of transfer function data. The program would be changed so 
that a denominator factor would not need to be entered by 
the user unless present in the transfer function being 
entered. If not entered it would always contain unity. 
This modification would be particularly useful when study- 
ing the effect of proportional gain on the performance of 
a control system. 

Another modification might be to include a step of 
unity as a standard forcing function. The system presently 
allows a unit step to the setpoint to be called for when 
entering the input variable of interest. 

The third possible modification would be to change 
the method of calculating the steady state value. The Final 
Value Theorem would be used at the time the overall system 
transfer function was determined instead of the present 
method of determining an exponential envelope that encloses 
the transient response curve. The method that has been 
implemented works quite satisfactorily as long as the re- 


sponse curve is quite oscillatory. 
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An addition that might be helpful would be the 
inclusion of the values of the poles and zeros of the 
open loop transfer function on the root locus diagram. 
These values are not determined by the present program 
but should not be difficult to include. Some of the 
other rules that appiy when constructing a root iocus 
diagram such as the asymptotic behavior, the angle of 
departure or approach to the real axiS, and the breakaway 
point could be implemented. 

The program could be modified to include the z- 
transform capability by modifying the transfer function 
entry section of the CSDAP system. Program logic avail- 
able in Lofkrantz's work (24) could then be used. A 
possible method of implementation would be to require the 
user to enter a z in front of the two transfer function 
name characters when it is desired to enter the transfer 
function as a polynomial in z notation. 

Another special type of transfer function entry 
similar to that for a time delay could be used for entry 
of zero, first, and second order samplers once the system 


was modified to handle z-transforms. 
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Appendix A 


USER'S MANUAL 


Ax lehelIntroduction 


This section is intended to be the user's manual 
for the Control Systems Design and Analysis Program 
(CSDAP). This manual is divided into five sections that 
describe the five major operating sections of the CSDAP 
system. 

The first is concerned with bringing the program 
into operation. After initialization it is necessary to 
enter configuration data for the control system to be 
studied and this is the concern of the second section. 
The last three sections outline the calculational pro- 
cedure for determining the transient response, frequency 
response, and root locus for the control system. 

Numerous tables and figures are provided through- 
out this description to clarify the specific operating 


procedures being discussed in each section. 


A.2.° Computer Operating Instructions 


The following description concerning the use of 
the CSDAP system will proceed assuming that all necessary 
programs and data files have been successfully loaded 
onto the system. If this is not the case, the Systems 
Manual should be consulted for procedures to follow in 


loading the CSDAP system. 
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A.2.1. ARuncInitiation 

This system as is described here is implemented 
on the IBM 1800 Data Acquisition and Control Computer. 
It is run as a Nonprocess Program under the IBM 1800 
Timesharing Executive System (TSX). To initiate the 
program it is only necessary to load the following 
cards into the card reader and depress the Console Inter- 
rupt button (same procedure for any other Nonprocess Pro- 


gram). 


4 10 15 20 40 

AL aeIO’ x x x USER NAME 
4 8 16 

// XEQ CSDAP FX 

LA wkSOB x xX X 
4 

// END OF ACU JOBS 


BLANK CARD 


A.2.2. Operating Instructions 

The next three subsections outline the operating 
procedures for the CSDAP system. They are general and 
apply to all phases of the CSDAP system and are given 
in this section rather than being repeated in the other 


sections. 
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A 272.246 Gontrol Digit Entry 


The control systems design and analysis program 
(CSDAP) is divided into nine logically self-contained 
sections. The entry into each of these sections is 
governed by the master control or executive program. 

This executive requests the user to enter the CONTROL 
DIGIT for the next section to be executed. A return 
is made to this executive after the completion of each 
section. Table A-1l lists the functions performed when 
each CONTROL DIGIT is entered. 

A listing of this table can be obtained at the 
user's typewriter by having data switch 1 in the "ON" 
position when the program is initiated into execution. 

At any time after that the listing of the CONTROL DIGIT 
functions can be obtained by having this same data switch in 
the "ON" position and entering a CONTROL DIGIT of 0. The list- 
ing may be terminated by resetting the data switch to the "OFF" 
position. 

A.2.2.2, Data Switch Functions 

The user may need comprehensive input help to pro- 
ceed to use this program the first few times. Later on 
this input help, depending on the sophistication of the 
user, May or may not be appreciated. If a user is very 
adept, excessive verbosity on the part of the program 


will be considered a severe annoyance and will make the 
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Operation slow. If the user is a novice or highly 
subject to making errors, the help is appropriate to 
insure proper input. 

The option to obtain assistance or not has been 
provided in CSDAP by means of data switches. The set 
of switches used depends on the setting of Sense 
Switch 2 on the IBM 1800 operator's console. If this 
switch is set in the "OFF" position, the console data 
Switches will be used but if it is set in the "ON" 
position it causes the set of digital input switches 
to be used. 

The data switches of the set chosen above when set 
in the "ON" position cause the program to list all 
possible alternatives or options available for the func- 
tion parameter before it is requested to be entered. 
Table A-2 lists the functions of each data switch. 

In general it can be said that anytime the user 
wants more than a single line of information when the 
input request is made, the pertinent data switch should 
be set to the "ON" position. 

All data switch settings are updated by the 
executive program immediately after a legal CONTROL DIGIT 
is entered. This means that there will be time for the 
user to correctly set these switches before proceeding 


to the next program section. 
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An example of the additional information that is 
typed out to the user is illustrated in Figure A-l. The 
first section shows the correspondence between the user 
and the program if data switch 3 was in the "ON" position. 
The same data entries are made in the second part of the 
figure but now data switch 3 is in the "OFF" position. 
This gives an example of the varying amount of support 


that is available to the user. 


A.2.2.3. User Data Input 


The program uses two subroutines written in assembly 
language for ease in entering the input data. The first 
of these allows entry of a single numeric digit from the 
typewriter keyboard. This was necessitated because many 
of the parameter options could be coded by using some con- 
tiguous string of numeric digits, 0-9. The end-of-field 
character does not have to be typed in addition to the 
desired digit therefore improving the speed of operation. 
A legality check for a numeric character within the desired 
range is made internal to the subroutine and input is 
requested again if in error. This is illustrated in the 
following statements taken from a typical run of CSDAP. 


The user's input to the program request is underlined. 
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4 = UNIT STEP -TO SETPOINT 
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ENTER INPUT VARIABLE FORCING FUNCTION 
(i Sa gaa 
2 = PULSE 
a5. = 9S INE 
4&4 = RAMP 
5 = JMPULSE 
2 
ENTER PULSE HEIGHT AND TIME DURATION 
Neo ee «0 
ENTER CONTROL DIGIT 
3 
PIGURE A=-1A EXTENSIVE USER ASSISTANCE FROM CSDAP SYSTEM 
ENTER CONTROL DIGIT 
2 
ENTER TRANSTENT RESPONSE OF INTEREST 
1 
ENTER INPUT VARIABLE 
1 
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Zz 
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FIGURE A-1B LIMITED USER ASSISTANCE FROM SCSOAP SYSTEM 
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ENTER TRANSIENT RESPONSE OF INTEREST 


0 = OPEN LOOP 
1 = CLOSED LOOP BEFORE FEEDBACK 
2 = CLOSED LOOP AFTER FEEDBACK 
ug 
INPUT IS IN ERROR, TRY AGAIN 
3 
INPUT IS IN ERROR, TRY AGAIN 
= 


ENTER INPUT VARIABLE 


The first error was caused because a non-numeric charac- 
ter was entered while the second error resulted because 
the numeric character entered was not in the range 0 to 2. 
The other subroutine allows the CSDAP user a less 
restrictive method of data input than would be possible if 
the program were written in FORTRAN. Integer and real 
variables in addition to alphanumeric data are read in as 
Single words or into vectors. The variables need only be 
separated by a space or comma instead of being right 
justified in their respective fields as required by a FORTRAN 
request. Leading and/or trailing blanks are ignored with 


the decimal point of floating point numbers being necessary 
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only when a fractional part is present. 

This subroutine also performs a check on the 
input data and results in notification of the calling 
program for suitable action,;usually to request the 
data again. The following few lines show the use of 


this subroutine with the user's reply underlined. 


ENTER INITIAL AND TERMINAL FREQUENCIES 


CSDAP 

, INPUT IN ERROR, TRY AGAIN 
re INPUT IN ERROR, TRY AGAIN 
-8 2.5 


The first error resulted because non-numeric data was 
entered for the data required. The last error was 
caused because only one of the two frequency values 


was entered correctly. 


A.2¢3%~, Run.Termination 


To terminate operation of the CSDAP system, the 
user simply enters a CONTROL DIGIT of 9 which causes 
the program to give control back to the Nonprocess 


Monitor of the IBM 1800 TSX system. 
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Avo. Transfer Function Specification 


Reo we (CONtrol System Representation 


The analysis of linear feedback control systems 
is facilitated by the use of transfer functions and 
block diagrams. The block diagram is a functional 
representation of an entire control system with the 
component transfer functions shown within the blocks 
and the flow of signals shown by lines connecting the 
blocks. Each block is an independent unit in that the 
connection between the blocks does not affect the 
transfer function inside them. 

A simple single loop control system is shown in 
Figure A-2. It serves to show some of the special trans- 
fer functions that can be calculated from the two blocks. 
Seldom is it possible to represent a control system with 
such a simplified block diagram. The control engineer is 
usually interested in analyzing the system performance for 
a number of input and output variables. Consequently, to 
handle the majority of systems that would be of interest, 
it was decided that the general system representation as 


shown in the block diagram of Figure A-3 would be employed. 


ee LOCK Diagram Transfer Functions 
To facilitate easy block diagram manipulation, the 


transfer function data for each block of Figure A-3 is 


OL=A | iv ; 5 


7a ya on ee 

Hobtsoiticoega aoxtai oka 39%: past Ged 

Y <) > ; wi a a 
noltezapgsiqed metaye Lozano” LEA 


7 


bos enotsonu? xe*ans213 to sam edt yd botssilios = 


sncicgonut 5s at msxpsidS Aoold oir * amsipelb ; 


‘2 
ofy diiw meseye Llougnes suisas 16 to noitsinee 


joie anoizoau® setens 23 tnenoamo 
‘eit partosnnos aeotl yd awode e ulsnpe. to wolt on 

: ‘ = x 
nmebnedebai ts ef Wold dost a 


sit tostts ston asob exncid id neswied mo itos * 


ke oo a 


- mond obite at nisdearasnet’ 793 eTs4 “ail 
~ Th — 7 : 

ait mwoda ei metaye.loxrtnoo gooh glence elqmia A * 
-ansit l[eiosge sit to smoe words oF Seuaah i > on Ser 


P s 


—— 
7 
, 


,exsold owz sot mort pods tos mice ed aS> ted3 enottani 
osu 4 
ijiw mateye Lloisgnos & Jnseverges és ofdseeog ob er 2 mob : 
- ' 
ei xsenipnas Loujnos siT .mpipSsth Aobd: Bordtignts | ‘a 
10% sonsmxoiieg meteye eft pass visrs wd besesterat y vitaue i 
* 


ot ,yfjneupsencD .asidasisv Sngaate bus, Yuqak hos 


; * ne 
,tesuestat to 5d Bieow tars am es8ve ro wtexotam ° 


es molssinesesitas: meteya f sxemae bape fe " 


C (s) 


NOMENCLATURE 
R(s) - Setpoint or reference input 
C(s) - Qutput or controlled variable 
G(s) - Forward path transfer function 
H(s) - Feedback path transfer function 

G(s) H(s) - Open loop transfer function 
G(s) 
1 + G(s) His) - Closed loop transfer function 
1 + G(s) H(s) = 0 - Characteristic equation 


FIGURE A-2 Simple Control Loop System 
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entered as a ratio of polynomial terms in the complex 
variable s except for a time delay. The method of 

entering a time delay is explained in a later section. 
A transfer function can be described by the following 


general expression. 


2 
(a.se+ az) (agSémbnaaSnt ac) deat) sa. 
ips) = See Ss eee se. To See (Red) 


2 
(b,s + bo) (bs + b,s + be) (.--) eee 


The transfer function is specified in the program 
by the two alphanumeric characters ID which can be any 
of those shown in Figure A-3 such as P2, Cl, or H2 to 


mention a few. 


A.3.3. Transfer Function Data Entry 


The coefficients for each polynomial of a transfer 
function, of the general form shown in Equation (A-1), 
are entered within brackets and appear ordered from the 
smallest to the largest power of s. The first value is 
always the constant term (coefficient of the zero power 
of s). The distinction between numerator and denominator 
polynomials is made by placing a slash or division sign 
after the closing bracket of the last numerator polynomial 
specification. 

Input data for specification of the transfer func- 


tion given as Equation (A-1) is entered by the user as: 
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ID (a5,8)) (€,,4,783) 22+/(b5rb,) (be ,by,b3) --- (A-2) 


To configurate a typical control system the trans- 
fer functions of a number of the blocks must be entered. 
The order of entering the block diagram transfer func- 
tions is irrelevant. It is important to realize that 
the program is written in such a manner that whenever a 
new transfer function is entered correctly it automatically 
replaces the earlier version of the same transfer function. 

The data is checked for a legal transfer function 
(block) name and the correct coefficient arrangement imme- 
diately after input. If the data passes the check it is 
transferred to permanent storage. If an error is found a 
Suitable error message is typed to the user depending on 
the type of error. 

If there is a name error the input section of the 
program will return to the executive to request another 
control digit. This provides the means for exiting from 
this section of the program. When an error in entering 
the input data occurs it is assumed that the user is seri- 
ous about entering this transfer function data and the pro- 
gram requests that the data be re-entered. 

Figure A-4 shows the information the user had to 


supply to CSDAP for configurating the control system 
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given below. 


where 
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For the entry of these transfer functions the user can 


let Cl(s) = G(s), P2(s} = G (8), and H1l(s) = H(s). 
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FIGURE A-4 CONFIGURATION OF A CONTROL SYSTEM 
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: ; -TS 
The pure time delay or transportation lag, e ‘ - 


is frequently found in control systems. It cannot be 
expressed exactly as a ratio of polynomials as was the 
case for the other transfer functions previously dis- 
cussed. It is therefore necessary to use an approxi- 
mation method. 

The program allows two approximation methods. One 


method involves the replacement of e °° 


by a rational 
function of s such that the first terms in the Taylor 
power series expansion of aig appear in the denominator 


as shown in Equation (A-3). 


i 
e = ee eee eee «6 (A- 3) 
e Lee a: Paso) + re / 3) eee 


In this program up to a fourth-degree Taylor series expan- 
sion can be used in the denominator. 

The other method is known as the Padé approximation 
method. In this case the power series expansion is approxi- 
mated by an equal number of poles and zeros. -.ne approximati 
improves as the number of poles and zeros used increases. 

In this program a Pade approximation with up to four poles 
and four zeros or in other words up to a fourth degree 
polynomial in both the numerator and denominator can be 


requested. The general fourth degree Padé approximation 
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is given in the following equation. 
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To enter a time delay approximation the user 
enters a transfer function name of TD when requested to 
enter transfer function data. The program then requests 
the user to enter further information describing the 
desired time delay approximation. 

The time delay can be placed in the forward path, 
as block P3, or in the feedback path as block Hl. Once 
the rational polynomial approximation is calculated it 
is multiplied by the existing contents of the respective 
block P3 or Hl. This makes it possible to have another 
transfer function already present for this block and dic- 
tates that the time delay should be the last transfer 
function entry for these two particular blocks. 

Two examples of the information that must be pro- 
vided by the user for obtaining a time delay approxima- 
tion in the control system are shown in Figures A-5A and 
A-5B. Figure A-5A shows how a third order Padé approxi- 
mation could be placed in the forward path while Figure 
A-5B shows a fourth degree Taylor series approximation 
being configurated for the feedback transfer function 
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A vaso. Data Input Checking 


Configurating the control system to be studied 
must be done with the most care. This is because the 
Operation of entering the transfer function data is 
most susceptible to errors. The input data is checked 
by the program for the correct transfer function name 
and the proper configuration of brackets. The coeffici- 
ents themselves cannot be checked because only the user 
knows these values. 

To facilitate additional checking of the parameters 
entered data switch 2,if set in the "ON" position, causes 
a transfer function to be immediately listed back to the 
user if it was successively entered. For the numerator the 
first line is an integer description that gives the number 
of polynomials that were found in the user supplied data 
plus the degree of each of these polynomials. The next 
line gives the coefficients for each polynomial successive- 
ly. The coefficients are ordered from the lowest to the 
highest power of s in a similar manner to the transfer 
function input data. The denominator data is listed in 
the same format. 

An example of this data listing is given in Figure 
A-6 which is the same transfer function input data as 
shown in Figure A-4 but because data switch 2 was set in 


the "ON" position, the additional listing of. the transfer 


i i Rat i 
} Leo 
Md ; _ Py ‘ 
‘ At, . 
Of ~A i fA 
I A 
-pattsodd, 3 uigal ets 2.t.A- 7 
cette nt ey a 
bute sd ot masays Loudnoo- ont enizexppr acd he,“ 


t 
ley 
2) 


25 6isb tugqnat aoisome’: setans13 wens en; 


7 


1er not 

tis09. S 
ot 4 

<A GS — 5 9 be § 


> f¢ 
a 


od +92. Bew S$ dositwe: Bt sb aaNead a0 b-s hes’ 


sb noisonui tstens13 sid privedne to aois 


dt .etextosid to aoldsig pitnen reqozq © “oa 


vino se2usDsd: bsAosno edd JORAES pres, 


tiaoq “WO” oft at tos ti 5 dot iwe. sisb ° 


simonyloq dose 104 a9 iestien tae eit 5) 


>t ath? -.srer Seon “ark faitiw oe 


- 


do et stab sugnit sdT vatorte of afdisqesseve tom 
: aa ‘ 


serrry ow As 
a \ 


gS 
Pe 2 


1s1t too1x0> sas 101 pate on: i < 


t 


_. aeuisy seott 


to prt wey edo Isno if: i tbbs ag Sgtlios? of 


= ~ ay : 
betetl y [636 tr banat - ed ot moeksong? sta: cand r 
‘ vi > is 


eit tot .bexstas ylevigessoue esw $4 +e 


. 


vip ted nobtgizoese ropes om Ss oaks send 


ssav sAdt at bavos s1sw tails elbimony fog 


: 


-eisimonylog seeds to rose: to ata heal edd 8 aut 


ol eft moxt bevebro 926 singioitives ® 


od3 o+ xomnem telimie 6 nte to. a 


© 


i sasb sotsnimoned of seb Juqak a 
- —- 7 -_ 
- } 


sivpit mi moevip et otcavat 6d Bh aids 20 ofa 


“SS rk 


- r 
_ _ 1s or. * 
re eee ae 
if a> 7, oie 
a es 


“~~ att Lat - 


A 
EN PER BGONFEROL «&hGI t 
a 
ENDER CONTROL (DIGHT 
1 
EN EER TRAWSE-ER-FUNGT LONe DATA 
Be a) tate) Se Lt 1 0s Soi Caty ts) 
2 0 ts 
22050000 1.10000 160000U 
3 1 1 : 
0290000 1.00000 1200000 0225000 1600000 
1.00000 


ENTER TRANSFER FUNCTION DATA 
ie co ils Oe 25.0 4N/ 9/00 
BRREORSIN@G TI NPUT@DATAL= ERYSAGALN 
te 2S) (160s 33040050033) 
2 0 1 
La2d000 1.009969 0.330C0 


is 
0200000 0233000 


Siler TRANSFER FUNCTION DATA 
eaer ) S190 % 54) 


1 e 
1200000 
iM 1 


1200000 Os 5UDUG 
ENTER TRANSFER FUNCTION DATA 
TELCEGAL- TRANSFER FUNCTION: NAME 


ENTER CONTROL DIGIT 


FIGIRE A-6 DATA VERIFICATION AFTER INPUT 


21 


function polynomial coefficients is provided. 


Peto Ore Listing of Transfer Functions 


If at some time the user wants to know the current 
status of some or all of the transfer functions in the 
system, perhaps because of modifications to the original 
control configuration, this can be obtained. 

This service is available by entering a CONTROL 
DIGIT of 7. The user specifies the output device, type- 
writer or line printer, where the transfer functions are 
to be listed and then enters the names of the transfer 
functions to be listed. If at any time the user specifies 
a transfer function name of 99, a complete listing of all 
the transfer functions will be made on the line printer. 
The transfer function names, and their relation to the 
black diagram of Figure A-3, are given in Table A-3. 

The data that can be listed on either output 
device is given as a single polynomial formed from com- 
ponent polynomials of the transfer function. Again the 
coefficients appear in ascending powers of the s variable 


with the 3? 


coefficient always appearing at the beginning. 
Figure A-7 shows the entries to list some of the 
transfer function data and the result while Figure A-8 


gives a complete listing of all the transfer function 


data. 
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A.3.7.: Initialization of the Transfer Function 


Area 


All transfer functions are initialized when CSDAP 
is brought into operation or later by entering a CONTROL 
DIGIT of 8. This procedure sets the degrees of the 
numerator and denominator polynomials equal to zero and 
sets the constant coefficients associated with zero degree 
polynomials equal to unity except in the case of the feed- 
back transfer function of the inner control loop, namely 
H2(s). Its numerator is set equal to zero to indicate 


zero feedback on this inner loop. 


A.4. Transient Response Calculation 


This section of the manual will describe three 
phases of the CSDAP that are utilized in performing the 
time response calculation. The first phase is used to 
calculate the overall control system transfer function 
and then to perform an inverse Laplace transformation of 
this transfer function. The second phase calculates and 
displays the time response data for the previously inverted 
overall transfer function. The third phase calculates and 
lists the performance criteria associated with the time 


response calculated in the second phase. 
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A.4.1. Transient Response of Interest 


The transient response of the variable of inter- 
est is the first information that is requested since the 
program allows the use of different transfer functions. 
The transient response of the open loop transfer func- 
tion, defined as the transfer function B(s)/R(s) with 
the feedback path on Figure A-3 broken at the summing 
junction, is one optional response that can be calculated. 
The input variable illustrated here as the setpoint R(s), 
can be the setpoint R(s) or either load variable. 

The response of the closed loop transfer function 
considered as C(s)/R(s) or B(s)/R(s), as defined in 
Figure A-3, can also be obtained. The option of speci- 
fying either closed loop transfer function is included 
because although it is common to refer to C(s) as the 
output variable, this variable usually cannot be measured. 
The control signal is that from the measuring element 
transfer function Hl(s), namely B(s). 

Figure A-l gives an example of the input to the 


CSDAP to specify the desired transfer function of interest. 


A.4.2. Input and Output Variables 
The open or closed loop transfer function is used 
to allow the design engineer to study the effects on the 


system behavior of changes in different input variables. 
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The procedure is to decide which of the possible 
inputs of the system to study and to call this the input 
variable. Another property will be called the dependent 
variable and may be considered the output. A relation- 
ship between the input and output variables is determined, 
and then the change in the input variable is introduced. 

The setpoint, load one, or load two are available 
as input variables while the output variable can be either 
C(s) or B(s), as discussion in Section 4.2. The different 
types of disturbance that can be used for any of the input 
variables are shown in Table A-4, which in addition to the 
codes gives the additional information required for each 
type of disturbance. Examples of the information supplied 


to the CSDAP are given in Figures A-3 and A-9. 


A.4.3. Inverse Laplace Transform Calculation 


When the user has specified the desired input 
variable and the type of forcing function with its addi- 
tional parameters, the program proceeds to calculate the 
overall transfer function. The inverse Laplace trans- 
formation is calculated from this transfer function by 
the method of partial fractions. The data vectors asso- 
ciated with the partial fraction coefficients can be 
optionally listed for the user by setting data switch 4 
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wa mn 
' 
i] if : 
I€-4 
r : j 
Le = t 
sidizsoq oft 33 doidw Shine od BL 


- 


Suani sit eiddt Liss of bois ybhuce od sie 


tnebriesaesb sit beliss sa. tiiw we reaqaa 1edtonk ved 6L=S\ 


-~noitsisx A .IvNGqsIHO ott Bsarsbianoo sed ysm bas eideias BV 


,benirmis355 ei esid&iitayJugiuo bas jugs odt preraieie. ide 
2 q4 a 
~bsoubousni et sidalasv tuqai sds icy spitsiy ent neds bas 


sidelisvs sits ows Bsol 16 ,SsH0. Beal iekogson ont 


ml 
jnexsttib ofAT .S.3 moitose at gobresvoetb as , {aja x0 ts 


7 


tugni sda to yos tot seu sd ass t6d3 soasdapselb to as 
re / p= in 

efj ot motitbbs at dotdw ,b-A sid6sT ai awone Se 
i. . 

fionss 101 bsxrleps1 noises otat isnots fhbs edt eov.ie 


betigqqne aottsmxotal sat To aolqnsexd -eonsdiutetb to ¢ 
+: } P 
“A bis €-A sesxzvpet re newitp 16 qAaaD 9 


ii 


motsealusotisd ozote 6 <T Sos tx Bi 985 eateV wad ee 


eit, eteluolso ot abssoorg ms2potg one -yahabel . 
; xaeik | ws 
-efexs soslaqsd sexrsvak srit arty ahoa enh 


yd notsoayvi xasteaness elds mozt Pavalvotes: poh: 


-O885 ertotosev 5i8b sdT - -AnOLs9 8x2 sateene 
ye 
, 


ed ms90 estneiottisos Hol ieest Ine: 


a tos iwe 6355 enketen ed teey ort Bet 
' Ah i a ‘*! 7 rs dd : 7 Le id A ¢ Wi: 

: is Mm ot & oe a) ‘et ££ 7 (Dp anh § hie & 
Prati el i A ey | 


LASLE 2A 


—— oe oe 


RAMP 


IMPULSE 


ai FORCING SEUNCT [ON 
Pikeot 
CODE DISTURBANCE 


NUMBER 


——<— a oe 


UI 


PARAMETER 


ee ne et ee ee ewe ee wes Oe ane Oe 


MAGNITUDE 


PULSE ane Gi) 


AMPLITUDE 


SLOPE 


CODES 


SECOND 
DISTURBANCE 
PARAMETER 


oo ot or oe Oe ee 


shal 


ENTER CONTROL DIGIT 
ENTER GCORTROL “OPGET 
ENTER TRANSFER FUNCTION DATA 
Ome) /ERSPOEIT.S) 


ENTER “TRANSFER FONCTION DATA 
Cie 0G) F023 ).(05))) 


ENTER TRANSFER FUNCT TOM DATA 
ACE / Sts te) ACT) 


ENTERSTRANSPER FUNCTION DATA 
TLLEGAL TRANSFER FUNCTION NAME 


EN? GR rCONFROUADEGRT 


a 
MVERSE TRANS RCRA TION FAS NOT BEEN CALCULATED 
SINCE «TRANSFER FUNCTION AREA CHANGED 
EGTER CONERG@L DIGIT 2 
ENTER CONTROL DOPGIT 
re 
ENTER TRANSIENT RESPONSE OF INTEREST 
ENTER INPUT VARTABLE 
i 
PNIER INPUT VARIABLE FORCING FUNCTION 
1 
ENTER STEP: MAGNITUDE 
1.0 
LINEAR ROOT CONTROL VECTOR 
2 A l 
LINEAR ROOT DATA 
A AQ Al A2 


—Gea loo ot OU, et video, “OsQUOUR. OO 
O.000CE CO O.1000E CO O.00GCTE OD 


QUADRATIC ROOT CONTROL VECTOR 
1 1 
QUADRATIC ROOT CATA 
A BR THETA 
—O0e4581E 00- O666459E=01 064643EF 01 


PIGURE SA] CALCULALION. OP SLHE INVERSE 
LAPLACE TRANSFORMATION 


UsVUCVE O0 
OCeOUDUVE O00 


C 
O-1002ZE O21 


THETA! 


74 ‘. : } 
ui 


| . - 
tote: seetD ical 


- *yuye weld 4 
I 
\p les 
ut tyaeceFT 2, Se 
fae by S90 De _ 


3 pyagvart 9otee ee 
Vienite yy (OGgr se 


seeAatt AIT ae 
(ri (2287 aes 


i 


5 5 f2aReahT RATA 


i 


» 
_ 
, 


oe, 
at °PagHAPT ASI IID ee 
jou Jonte> saya. ae 
(i oo ae, Sepae an 
“8 MWeaPt S20te ae 

i re l ~ ie eva? egTi2) ) = 
s. ; - 

‘(gio seeTpes PaTeae 


— SS 
‘= 24eR4erT S3Goue 
Pe caw TUSK) OATRae 
‘shi PAS, TU BOTs 
& ‘ ’ 
scutieda te # ; 
thes ae 31V9 4 


4 Rae m, 
ceTug? Tog sean 


A-33 


the user is informed of the number of roots that exist, 
their type and multiplicity. 

An example of the output is given in Figure 
A-9. The linear root control vector informs the user 
the number of linear roots and the number of times each 
is repeated. The coefficients associated with Equation 


(A-5) are given under'the linear root data. 


F(t) =e (AO + Al*t + A2*t*t) + H(t) (A-5) 
where s = A is a linear root repeated n-1l times 


and F(s) = ee ase re das WG ig ie cael Rear 

(S-A)” 
The presence of quadratic rootsresult in a similar dis- 
play of data. If the quadratic factor appears only once 
in the overall transfer function it takes the form of 
Equation (A-6) while it assumes the form of Equation (A-7) 


if it appears repeated once. 
F(t) =Ce SIN; (B¥te+- THETA) + H(t) (A-6) 
where s = A + iB is a single quadratic factor 


and: F{s) = tis) cps he® 
(S-A)"-+ B 
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F(t) =e {C(SIN(B*t + THETA) -B*t*COS (B*t + THETA) ) 
+ C'*COS (B*¥t + THETA') } (A-7) 


where S = A + iB is a repeated quadratic factor 


and F(s) = i), + H(s) 
b(SsA)2" 2-RoBS) 


A.4.4. Transient Response Data Calculation 


The partial fraction coefficients must have been 
previously calculated before the transient response can 
be determined. If any modification to the block diagram 
transfer functions has taken place since the last par- 
tial fraction coefficients were calculated, an error will 
result when this phase of the program is executed. This 
error message is shown below and also appears in context in 


Figure A-9. 


INVERSE TRANSFORMATION HAS NOT BEEN CALCULATED 

SINCE TRANSFER FUNCTION AREA CHANGED 

ENTER CONTROL DIGIT 2 

Once the inverse Laplace transformation with asso- 
ciated partial fraction coefficients is determined, the 
program then requests the user to supply parameters about 


the transient response calculation. The transient res- 


ponse replot digit, the total time, and the number of 
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data values to be calculated are included as the para- 
meters. The last two parameters are only required if 
an Original plot of the transient response data is re- 
quested as given by the replot digit. (Figure A-10) 

A maximum of five hundred transient response 
data points can be calculated. The time axis always 
starts at zero time and proceeds to the total time. The 
response data is stored into a disk file immediately 
after it is calculated. Because of this feature the 
user can display this transient response data at a later 
time, by means of the replot option, without recalcula- 
eLOn. 

Figure A-14 gives some of the errors that can 
occur when entering the parameters necessary for the 


transient response data calculation and display. 


A.4.5. Display of Transient Response Data 


After the transient response data has been cal- 
culated and stored on disk, information concerning the 
display of the data is requested. This input includes 
the display device to be used, scope or plotter, the 
plot direction, and the title for the display. Four 
possible graph types, as listed below and shown in 


Figures A-1ll to A-13, are available. 
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i) Vertical Plot - Figure A-1l 


ii) Horizontal Plot - Figure A-12 


iii) Top Half Page Plot Figure A-13 


iv) Bottom Half Page Plot Figure A-13 


The transient response data can also be listed 
on the lineprinter. This is accomplished by having data 
switch 6 set in the "ON" position when this phase of the 


program is entered, (Table A-5) 


A.4.6. Performance Criteria 


A number of performance criteria can be calculated 
from the transient response such as rise time, steady state 
value, and percent overshoot to mention a few. These per- 
formance criteria are closely associated with the transient 
response of the system and can therefore be displayed with 
it. If a display of the criteria is requested it appears 
on the bottom half of a plot of the type shown in Figure 
A-13 if the results are displayed on the scope or the 
plotter. Another option allows the user to get a listing 
of the criteria on his typewriter. 

Figure A-14 shows the calculated values for the 
different performance criteria while Figure A-15 shows this 


same information as well as a plot of the transient response. 
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FIGURE A-13 Example of Two Plots Per Graph 
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AOS . Frequency Response Methods 
ATH5e15 Frequency Response Options 


Frequency response methods are based on the open- 
loop transfer function G(s) H(s) as given in Figure A-l. 
Plots are made by replacing s by jw in the transfer func- 
t20n, G(s) Hs). 

Three methods of plotting this function are most 
useful in the frequency response analysis and design of: 
feedback control systems and will be illustrated in the 
following sections. 

Upon entering the frequency response analysis 
section of the CSDAP system, some general parameters must 
be entered ky the user. These include the frequency re- 
sponse replot digit and the type of frequency response 
diagram desired, that is Bode, Nyquist, or Log-Modulus. 
If a replot of data existing on disk has been requested, 
the CSDAP system proceeds directly to the phase pertain- 
ing to the type of plot desired. If a replot was not 
requested, the user is asked to enter the initial and 
terminal frequencies for the range of interest and the 
number of points desired. The points are distributed 
evenly over the complete frequency range. The program 
then proceeds to the correct phase to display the fre- 


quency response diagram desired. 
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A.5.2. Bode Plot 


The Bode plot, which is also known as a Corner 
plot, is composed of two separate graphs. The first 
graph is a semilog plot of the magnitude, in decibels, 
versus the frequency, in radians per second, on a 
logarithmic scale. The second graph is also a semilog 
plot with the phase angle, in degrees, being displayed 
versus the frequency on a logarithmic scale. 

The user needs only to specify the device to be 
used to display the plot and the desired title for the 
plot in order to obtain a Bode plot. Figure A-16 shows 
typical user supplied information while Figure A-17 is 
an example of a Bode plot. 

The gain margin, gain crossover frequency, phase 
Margin, and phase crossover frequency values are given 


in the lower left hand corner of the two plots. 


Ar 5egds Nyquist Plot 


The Nyquist plot sometimes known as the Polar 
Plot displays the magnitude in decibels versus the phase 
angle in degrees on rectangular coordinates as the fre- 
quency is varied from the initial to the terminal value 
as specified by the user. 

The user can specify the scaling of the plot or 


use automatic scaling. The plot device must be selected 
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and if desired a title for the Nyquist plot to be shown. 
Specifying automatic scaling ensures that the axes will 

be labelled large enough to accommodate all data values. 
Figure A-18 is an example of a Nyquist plot with auto- 
matic scaling. A user specified scale allows the display 
of a small section of the plot in more detail such as 

shown in Figure A-19. The user supplied information to 
obtain a Nyquist plot is shown in Figure A-16 for both types 


Of scaling. 


A.5.4.  Log-Modulus Plot 


The log-modulus or gain phase plot displays the mag- 
nitude in decibels versus the phase angle of G(jw) H(jw) 
in degrees in rectangular coordinates with frequency as 
the parameteron the curve. The plot is a continuous curve 
plot with each value shown by a small plus sign as can be 
seen in Figure A-20. 

In addition to the log-modulus curve, some frequency 
response specifications are also listed at the bottom of 
the graph. These include the peak resonance and resonant 
frequency in addition to the gain and phase margins and 
phase crossover points that can also be displayed with the 
Bode plot. 

As can be seen from Figure A-16, for the log-modulus 
plot, the user selects the plot device and then specifies 


the appropriate title. 
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FIGURE A-17 Example of a Bode Plot 
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A.5S U5. “Tabulation: of Frequency Response Data 


A listing of the frequency response data on the 
line printer can be obtained whenever the frequency 
response section of the CSDAP is entered. This optional 
listing is controlled by data switch 10 which is set to 
the "ON" position if a listing is desired.(Table A-6) 

The listing option is independent of the particular 
frequency response plot desired and the replotvoption Sy sLt 
a Bode plot is desired, then the program uses the same 
title as for the Bode diagram. If the requested plot of 
the frequency response data was not a Bode diagram, then 


the program will request an appropriate title. (Table A-7) 


A.6.. Root’ Locus Analysis 


The wot locus data is displayed on a plot device 
aS a point plot. Both positive and negative values of the 
gain are shown on this plot. They are differentiated by 
plotting frequency points with positive gains as small 
plus signs (+) while those points with negative gains are 
plotted as small crosses (x). 

The user must specify the replot option digit (same 
as for time response or frequency response), the initial 
and terminal frequency values, the plot device and the 


title of the plot (Figure A-22) as shown in Figure A-21. 
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A listing of the data values including the fre- 
quency coordinates and gain values can be obtained on 
the line printer as shown in Table A-8 by having data 


Switch 14 set in the "ON" position. 
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A.7. Program Limitations 


In the design of the CSDAP system a number of re- 
strictions were imposed on the program due to core size 
limitations and arbitrarily chosen criteria. 

Transfer function data can be entered for any 
block given in Figure A-3. The symbols used for this 
data entry must be the same as those given within the 
blocks of Figure A-3. The CSDAP system can contain a 
Maximum of two hundred coefficients including both those 
of the numerator and denominator of all the transfer func- 
tions listed in Table A-3. 

In the calculation of the transient response of 
the control system, only arbitrary disturbances such as 
a step, pulse, sine wave, ramp, or impulse may be used as 
the forcing function. The characteristic equation denomi- 
nator is restricted to containing no more than twelve linear 
roots and eight quadratic roots. Each linear root may be 
repeated twice while a quadratic root may appear in the sys- 
tem repeated only once (two equal roots). When the above 
rule is violated a suitable error message will be outputted 
to the user indicating that a linear or quadratic root nea 
repeated n times. It is usually sufficient to solve this 
problem to reduce the order of the Padé approximation. 

The plotting programs in CSDAP only allow one set 
of calculated values to be plotted on a single graph. 

This restriction of no multiple plots is viewed as a 


minor inconvenience because the CSDAP was designed 
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as an interactive system giving the user a quick view of 
the resultant plot. Enough other numerical data is 
given for a more exact comparison of various results. 

The Bode diagrams produced by CSDAP are restric- 
ted to a phase angle in the angle between -720 and +720 
degrees. If the phase angle should lie outside this 
range, a discontinuity of one cycle of 360 degrees will 
appear on the diagram. 

The CSDAP system does not allow execution to be 
continued after it is stopped. Each time CSDAP is exe- 
cuted it initializes the transfer function area and all 
other system parameters. Replots of transient, frequency, 
and root locus response data that is stored on the disk 
may be obtained without re-entering the transfer func- 


tion data. 
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